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Es werden die Matrixelemente fiir Absorption bzw. Emission von Pho- 
nonen berechnet, wobei die Elektronen nicht durch die Blochzustande beschrie- 
ben werden, welche dem undeformierten Gitter entsprechen, sondern durch 
Wellenfunktionen, die in geeigneter Weise den durch die Gitterschwingungen 
bewegten Ionen folgen, derart, daB sie in deren Umgebung praktisch ihre dem 
ungestörten Gitter entsprechenden Werte beibehalten. Zur Beschreibung der 
Kopplung wird der Nordheimsche Ansatz und der Deformationstensor des 
Gitters verwendet sowie die Rückwirkung der Elektronen berücksichtigt. 
Für Normalprozesse stimmt das so gewonnene Resultat mit dem Bardeens 
(1937) überein bis auf Terme, welche infolge einer schärferen Approximation 
des Bardeenschen Matrixelementes auftreten und eine schwache Ankopp- 
lung der transversalen Gitterschwingungen beschreiben. 


Bekanntlich beruht sowohl der Ansatz von Bloch und Bethe für das 
Potential der Wechselwirkung eines Elektrons mit einem gestörten Kristall- 
gitter (Ansatz der ,,deformierten Ionen‘) als auch der Nordheimsche Aus- 
druck für diese Wechselwirkung (Ansatz der ‚starren Ionen‘) auf einer Ent- 
wicklung nach der Verschiebung der Ionen aus der Gleichgewichtslage, welche 
man nach dem linearen Glied abbricht mit der für thermische Gitterdefor- 
mationen sicherlich unproblematischen Begründung, daß die Elongationen der 
Ionen klein sind im Verhältnis zu der Dimension einer Wigner-Seitz-Zelle!). 
Man berechnet so Übergangswahrscheinlichkeiten zwischen den Elektronen- 
wellenfunktionen des ungestörten Kristalls und erhält in der Näherung quasi- 
freier Elektronen mit kugelförmigen Energieflächen das bekannte Resultat, 
daß die transversalen Gitterschwingungen nur über Umklapp-Prozesse an die 
Elektronen gekoppelt’ sind. Zu diesem Ergebnis kommen auch die verbesserten 
Methoden, welche den Einfluß der Gitterdeformation auf die Elektronen- 
verteilung berücksichtigen ?)?). 
Wie Blount®) betont hat, erheben sich zunächst Einwände gegen eine 
lineare Entwicklung nach den Verschiebungsvektoren der Ionen, wenn, wie 


1) Zusammenfassende Darstellung bei J.M.Ziman, Electrons and Phonons, Oxford 
Clarendon Press 1960. 
2) J. Bardeen, Physic. Rev. 52, 688 (1937). 


3) M. Bailyn, Physic. Rev. 117, 974 (1960). x Die 

4) E. J. Blount, Physic. Rev. 114, 418 (1959). EC is 
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etwa im Fall einer dem Gitter aufgeprägten Ultraschallwelle, die Kleinheit 
dieser Verschiebungen nicht mehr gewährleistet ist. Unter diesen Umständen 
bietet sich die Verwendung des Deformationspotentials an, welches nicht nach 
den Verschiebungen, sondern linear nach den Komponenten des Deformations- 
tensors €,, entwickelt ist. Die Anwendbarkeit des Deformationspotentials 
erscheint durch den Beweis der Äquivalenz zu den direkten Methoden für 
Halbleiter besser fundiert) als für Metalle, bei welchen die Vernachlässigung 
von Gliedern der Größenordnung k, k, &,. nicht unmittelbar zu rechtfertigen 
ist®). Andererseits erfordert der von Peterson und Nordheim’) vorgebrachte 
Einwand gegen die Berechnung von Übergangswahrscheinlichkeiten zwischen 
den Bloch-Funktionen des ungestörten Gitters Beachtung, besonders im 
Hinblick auf große Verschiebungen der Ionen, welche ja, wie man auf Grund 
der adiabatischen Näherung!) erwartet, die Elektronenwellenfunktionen in 
nullter Näherung ,,mitnehmen“, so daß man eigentlich Übergänge zwischei 
diesen ‚„‚mitgenommenen‘“ Funktionen berechnen sollte. 

Im folgenden wird versucht, diesen beiden Einwänden nachzugehen und 
das Matrixelement für die Elektron-Gitter-Wechselwirkung auf Grund einer 
Entwicklung der Störung nach &,. für Übergänge zwischen den im Sinne der 
adiabatischen Näherung ‚‚mitgenommenen“ Wellenfunktionen berechnet, wo- 
bei die Beeinflussung der Elektronenverteilung durch die Gitterdeformation 
durch ein geeignet modifiziertes selbstkonsistentes Bardeensches Verfahren?) 
berücksichtigt wird. Dabei ergibt sich in linearer Näherung in ¢,, das bekannte 
Resultat?) für Normalprozesse und eine Kopplung an die transversalen Gitter- 
schwingungen, welche jedoch für freie Elektronen (m* = m) verschwindet und 
implizit bereits in den früheren Rechnungen!) vorhanden ist, dort aber ver- 
nachlässigt wurde*). 


I. Formulierung und Transformation der Schrödinger-Gleichung 
Für ein Elektron im ungestörten periodischen Potential lautet die Schrö- 


dinger-Gleichung 
oy j 
— + Volt) y = (1) 


Die Eigenzustände yr = ur (r) e'!" sollen im Grundgebiet von N Zellen des 
Volumens 2, normiert sein, die zugehörigen Energien im Leitungsband seien 
durch E,(f) bezeichnet®) im reduzierten Zonenschema. Für das periodische 
Potential V,(rt) nehmen wir folgende Darstellung an: 


Vo(t) = & — 9). 
Dabei durchläuft g alle Gitterpunkte, und es ist 74% r 
fo — g) + Ui(t— 9), wenn r in Z,, 
|0, wenn t außerhalb Z,, 


5) J. Bardeen u. W. Shockley, Phys. Rev. 80, 72 (1950). “ a 
6) S.C. Hunter u. F. R. N. Nabarro, Proc. Roy. Soc. A 220, 542 (1954). 
7) E. 5 Peterson u. L.W. Nordheim, Phys. Rev. 51, 355 (1937). 
*) Anmerkung bei der Korrektur: Inzwischen wurden auf der Prager Halbleiterkon- 
ferenz Sept. 1960 ähnliche Betrachtungen von G. D. Waitfield vorgetragen [G. D. Whit- 
field, Phys. Rev. Lett. 5, 528 (1960)], Phys. Rev. 121, 720 (1961) 
8) Es braucht hier nur das Leitungsband betrachtet zu werden. 
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wobei Z, die um den Gitterpunkt g gelegene Wigner-Seitz-Zelle bezeichnet De & 
und v(t — g) das reine Ionenpotential, während U,(t—g) der Anteil dr 
Leitungselektronen am Wigner-Seitz-Potential im Sinne der Hartree- 
Fock -Näherung ist. BT 
Mit?) setzen wir 


| 
U, (t-9)=5- für |c—gl|<r,, 


wobei r, der Radius der sphärisch approximierten Wigner-Seitz-Zelle ist. — 
Nun sei die Periodizität des Gitters durch eine Schwingung der Frequenz A 
w gestört, so daß ein Ion, dessen Ruhelage der Punkt g ist, in den Punkt 


1 h 
a ER i(wg— wl) -i(mg— wi) 

verschoben wird. Dabei ist e der Polarisationsvektor der Gitterschwingung. 
und M die Ionenmasse, by, ist ein Amplitudenfaktor. 

Im so gestörten Gitter schreiben wir die Schrödinger-Gleichung in der - 
Form 


MAP a, (t)) Ve(t,t) = 

4 
wobei V,(t,t) ein noch unbekanntes Störpotential ist Nun transformieren 
wir diese Gleichung ein Koordinatensystem, welches sich mit der Schallge- 


wre 


schwindigkeit c = mis tw gegen das Laborsystem bewegt. Diese Transforma- _ 
tion lautet 


Dann finden wir 


- t) +2 Pt) + - PE, 


Nun müssen wir darin noch R,(t) in den neuen Koordinaten ausdrücken. 
Offenbar ist 


=9-ct+ 
=9+ 37, 


so daß wir für (7) schreiben können 


8 


Nun greifen wir die von Blount®) vorgeschlagene Transformation auf krumm- 
linige Koordinaten auf, indem wir zunächst ein kontinuierliches Verschiebungs- 
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Daraus folgt in linearer Näherung in : 

a wobei = 


2 Nun transformieren wir (11) auf die neuen Variablen r’ durch 


feld $(r’) definieren: 


st) = elb,, ett b* e-iwr) 
vw + b% 


Es gilt also speziell 


Nun gehen wir zu den neuen Koordinaten 


über. Zur Umrechnung der Schrödinger-Gleichung auf diese Koordinaten 
brauchen wir den zu T’ gehörigen metrischen Tensor g,.. Nach (10) ist 


o 


7 


Die Umrechnung der Schrödinger-Gleichung (8) auf die Koordinaten 7?’ 


ös 


liefert bis auf höhere Terme in — und rv | ; oe 
(11) 
—g + 8(t') — +V,(t') =ih a” 


= 


erhalten wir 
+ 8(r’ — ct) — 84 (t) —g) + 
a 
wobei 
1 / h = 
8(r—ct) = | e {Dy 4 (13) 
2Mo 
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Da nun U(r — 9) jeweils nur für rally —g| <r, von Null verschieden ist, onl ve R 
man, wenn die Wellenlänge 2 = “” hinreichend groß gegen den Zellenradius r, 


ist, mit (2) näherungsweise 


(15) 


linaten 


ist —cl)e 
eo 


+ — cl) — (0)} grad + = 


Sie ist das Ergebnis einer Reihe von Transformationen der a 
Variablen in der ursprünglichen Schrödinger-Gleichung (6). Läßt man in zn) 
(15) überall den Strich an den Variablen fort, d. h. identifiziert man sie mit 
den in der Ausgangsgleichung (6) stehenden, so erhält man (15) in anderer _ 
Deutung als Ergebnis einer unitären Transformation e‘®, welche den Hamil- | 
ton-Operator 


H=- +2 U(t — g — + Volt, 


= 


ten 


aus (6) in die der Gleichung ae entsprechende Form 
+ A(t) an 


és 
gebracht hat, wobei man e'* so = hat, daB H’ linear in = fund 


fm: ist. Für Vj liest man aus (15) ab: zo 


lt 


+ — — SW} grad 9). 


Wir legen den folgenden Betrachtungen diese letztere Deutung zugrunde, 
indem wir Vj + V3 als Störung auffassen werden, welche zu Übergängen 
zwischen den Blochfunktionen yr(t) = ur(t) e'!" Anlaß gibt, die die Eigen- 
funktionen nullter Näherung des durch (15) beschriebenen Problems sind. Das 
bedeutet, daß die für eine Störungsrechnung geeigneten Ausgangsfunktionen 
des ursprünglichen Problems (6) nicht die Blochfunktionen, sondern die 
Funktionen 


[7 
— — — 
Th Xa (t 


o 


6x, Ox, 


= *S we = yr — 8 (ut (19) 
(13) sind, wobei mit u(r, £) diejenige Funktion gemeint ist, welche auf Grund von 


r=r —3(r’—cl?) 


t'—ct als Funktion von r ausdrückt: 


(13) 


(14) 
_ | 
4 
B 
5 
Fe 
| 
> 
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a Die Funktionen (19) schmiegen sich also in gewisser Weise der Gitterdeform. 
tion an, indem sie mit dieser ,,mitgehen‘‘, und zwar derart, daß sie in der Nähe 

der verschobenen Ionen im wesentlichen ihre alten Werte beibehalten, analog 

der von Peterson und Nordheim’) diskutierten Vorstellung. 


II. Berücksichtigung der Deformation der Elektronenverteilung 


Zur Bestimmung des unbekannten Störungspotentials V5(r’, £) verwenden 
wir ein Verfahren, welches sich eng an die Bardeensche Methode?) anschließt, 
allerdings ebenfalls ohne Korrelations- und Austauscheffekte?) zu berück. 
sichtigen. 

Wegen der aus (2) und (15) ersichtlichen speziellen Form der Zeitabhängig. 
keit der BE H,(t) erhalten wir lege artis für die Ubergangswahrschein. 


lichkeiten W}+ zwischen den Zuständen (19) den Ausdruck 
Wii = "2 E(t’) + ho), (21 


so daß wir uns bei der störungstheoretischen Behandlung der Schrédinger. 
Gleichung (15) auf eine statische Störung H, = = H,(0) beschränken können. 
Die exakte Lösung von (15) entwickeln wir nach on Eigenzuständen (19) 
; und erhalten in erster Schrödingerscher Näherung 


n’ |H,(0)| En) 


Hier bedeutet n den Bandindex, den wir jedoch später wieder fortlassen 
können. 
Wegen der zu erwartenden Auswahlregel 


wobei = 22 § und § Vektoren des reziproken Gitters sind, können wir auch 
schreiben 


Pen = Yen +23 { BE 9, n’ Yt+ w+ Rn’ + Br gn’ Yt-w-8,n'}- (23) 


Daraus ergibt sich die zugehörige Ladungsdichte der Leitungselektronen 


o’(t), wenn wir 


(t) va, H+ ine + Bree te e]] 
ar Ist nun ®’(r) das von dieser Ladungsverteilung erzeugte Potential, und ist 


= y(t), (25 


(24) 


UF n(t ) Utn(t 


setzen und linear in Big rechnen zu 


3e? 


ar, 


+ 8(t) —3,—g]’ 


e(r + 8(t) — 8, — — 9)? 


I>r 


ic '» 
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wobei Dar 
|_| 
|_| 
Pq(t) = ( 
Gt 


forma. 
r Nähe 
analog 


enden 
hließt, 
‚erück- 


längig- 


schein. 


(21) 


nger- 
Önnen. 
n (19) 


(22) 


lassen 


a H. Stolz: Zur Theorie der ae Te Te in Metallen 359 


so berechnet sich V3(r) aus 
Vale) = (r) — Bike) 


Nun können wir schreiben 


= g(r — g) + yf (x), 


e(t— 9) 


&(t — g) 


($(t) — $) für !r—g] < ry, 


Da nun, nach der Grundannahme der ganzen Methode, 
A (rx —g) =— yt (t) yr (t) 


gelten muB, berechnet sich V3(r) aus 


47 
(30) 
Wir finden od 
Aus (22) und (23) folgt nun a “ 


Brg = — 


<f+ w+ &| Vi + RAR, 
E,(t + to) — 


Führen wir noch die Bezeichung 


1 ine Big+B 


Daraus folgt 


{Bia + Bra 721 —E,(t+w) 


By 


ine 1 1 


> 


wobei 
= 
und 
_ 
| 2 
(23) 
(24) ein, so finden wir mit (31) unter Beachtung von (24) . 
| 
nd ist 
(25) 
_ 
(26) 
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mit Bardeens?) Funktion 
1 1 Ly 


Gere Den Ausdruck (32) hätte man nun in (31) einzusetzen und alsdann 
+ w + & | V5| £> zu berechnen. Vorher ist es jedoch notwendig, die in (32) ein- 
gehenden Matrixelemente {+ w+ + V3|D zu untersuchen und, 
soweit möglich, die verbliebene Summation über f auszuführen. 4 


III. Berechnung der Matrixelemente etal 


Bei der Berechnung von <f + w+ &|V1|f nach (18) treten zunächst 
Integrale vom Typ 


eco(f+ w+ R, oz, t ) 


auf. Um dieses und weitere Integrale möglichst weit ausrechnen und auf in 
der Leitfähigkeitstheorie gebräuchliche Ausdrücke reduzieren zu können, ist 
es zweckmäßig, sich zuerst auf Normalprozesse (8 = 0) zu beschränken. Dann 
hat man 


2 
© Us 


IX, (f+ w, f) = )w(r)d’r+ öx,öx, Br 
(35) 
h 


Obwohl wir im vorigen Abel u an spezieller Stelle (24) die Funk- 
tionen yr(t) durch ebene Wellen angenähert haben, wollen wir hier mit mög- 
lichst geringen Forderungen an die Wellenfunktionen y;(t) auskommen. We ir 
verwenden daher die Näherung der effektiven Masse®) und setzen 
Uf w(t) = uf (t) + w - grade uf (r) 
Mit f = £+ w erhalten wir 


Bei Kristallen mit NEAR gilt 


ur (t) = u_r(t) = 1). 


| ub (x) 


Also ist 
* 2 f 
grade ur ax, éx, 


Nun ist (36) das Matrixelement des hermiteschen Operators ;, 4, , geht also 


in sich über, wenn man f mit f’ vertauscht und zum konjugiert Komplexen 


®) S. J. Pekar, Untersuchungen über die Elektronentheorie der Kristalle, Berlin 1954, 
Akademie-Verlag. 
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übergeht. Daher muß der in f' — f lineare Term verschwinden, die höheren 
Terme berücksichtigen wir nicht. 
In ähnlicher Weise zeigt man, daß 


- 
Nun ist 


| uf (t) -=— | uelt) (88) 
aber er ET 27 


so daß in Entwicklung von nach f’ — f das in f’ — f lineare 
Glied verschwinden mus. Zusammenfassend erhält man also für (35) 
+E fut (i # 
cE, _ 


f(r) =) k oka + Des) — k, 


und der bekannten Relation 
J Ur u t= he (40) 
Unter Ausnutzung der kubischen Symmetrie man leicht (s. Anhang), 
= 
C( (f) (41) 
ist mit ‘ea Betheschen Abkiirzung 
h2 
(uf. + ur Auf ) dr = C (f). (42) 


Da nun nach (2) 


= b,eiwr _ H*¥ o-iwr = 
eo(t) VN ; Mw 2 Ow 
und 
h 


erhält man mit (39), (40) und (41) zunächst 


J2 (e, (e f) (fw) + (ef (m grad: Ey) (43) 
VN 2Mo|3 

bi] 


(33) 
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Als nächstes hat man zu berechnen en ti 
Jf — 30) grad U(t —g) yt yr 
8 


= fee) grad Vo(t) yf yr — grad Vo(t) pit yr dr. 
PR 


Hierin ist das erste Integral das Matrixelement der Wechselwirkung nach 
Bloch und Bethe, daszweite das entsprechende Matrixelement nach J ones”). 
Nach einer neueren Rechnung von H. Puff!!) erhält man mit (9) bei Be- 
schrankung auf Normalprozesse und in linearer Näherung in tw mit der Ab. 


kürzung re grad; E,(f): 

a fae) grad Vo(t) yp dr = \* 6a," 
+ [Bo(t) — + w +> + (eB [0,2 (ww 


Beachtet man nun (41), so erkennt man, daB 


Div 
+ (e (xt) +(e x) (Ew) + (Ey — Ey(f + mw))- 


Unter Verwendung der Schrödinger-Gleichung zeigt man '), daß für Normal- 


prozesse 
b 


—i (te) + w)) 


Nun ist 
i . 
+ if y(t). 


Also bleibt: 
8 a 


b 


10) H. Jones, Theory of Electrical and Thermal Conductivity in Metals. Handbuch 
d. Physik (S. Flügge), Springer- Verlag 1956, Bd. XIX. 
1) H. Puff, Ann. Physik 7, 368 (1961). 
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v 
Ve 
D 
V 
fo 
(44) 
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| fi 
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u 
fi 
F 
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und linear in mw: 


3, [grad V, By Va -{ 
2 je (t) pr yr VN 


iJ be] 


+ [Bg(t) — + w)] (e, ©) — 
es 10) (x) = G(wr,) = 3° cos sinz 
r Ab- Vereinigt man nun (43), (44), (45), so findet man _ eee Sa am a 
i f vere 
= — —E G 
+ | Volts) — Ho(0)] - 
Dw (47) 
\ 
J Dieses Resultat erlaubt folgende SchluBfolgerungen: LaBt man dic durch 


V5 beschriebene Wirkung der Deformation der Elektronenverteilung zunächst 
fort, so unterscheidet sich das Matrixelement (47) vom gelaufigen!) durch 
Terme, welche, wenn man den Energiesatz 


(44) Et) — Ey(f + w) = 
der in die Ubergangswahrscheinlichkeit eingeht, beachtet, in der Grenze freier 


Elektronen verschwinden, wenn namlich a =f wird. Für nicht völlig 


mal- 

freie Elektronen sind diese Terme von der Größenordnung h w (1 _ 2) (e, f). 
| Ihr Auftreten in der Übergangswahrscheinlichkeit ist keine Konsequenz der 
s unitären Transformation (17). Rechnet man sie unter Verwendung von (45) 


für Übergänge zwischen Bloch-Zuständen aus, so kommt man zu demselben 
Resultat wie mit (47) für die Übergänge zwischen den Zuständen (19). Wo 
man jedoch die Matrixelemente (45) beziehungsweise (47) ohne ö(E,(f + tw) — 
E,(f) X h w) braucht, wie etwa im Bardeenschen Verfahren der Berechnung 
von V3, wird ihre Verschiedenheit wirksam und sind Differenzen in den Resul- 
taten zu erwarten. ev 

Wenden wir uns nun der Berechnung des Matrixelementes 


zu. Beachten wir (29) und (28), so sehen wir, daß 


> [grad [Uc(e — g) + Us(e — g)] {8(t) — 8(g)} ddr 
) k 


* 


Aus ?) können wir sofort entnehmen: 
af 
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[Vo(r,) — G(x) 

5 


= > eit-Pi ma [grad et fe grad = 0 
q 
2 2 


für Normalprozesse in der Näherung der effektiven Masse, wenn man die 

Kugelsymmetrie von U, beachtet. Ebenso sieht man, daß das Integral ent- 
va . 

sprechender Art fiir U? verschwindet. Also bleibt: 

i f ine 1, iJ bu] 

(e, mw) -—G(zx) (49) 


YN 2Mw w 2 1-58] 


Für das folgende ist es zweckmäßig, noch (47) und (49) zusammenzufassen: 


VN '2Mo 


[em G(2) [Vor | 


Q, w? 


Geht man damit in (32) ein, so ergibt sich 


| h by 4 


N 
Ey(t) — y(t + w) By(t) — — w) 


1 ST _ e in 
(eb mold B,(t) — + w) 2 


f 
Da nun grad; E,(f) grad, E,(— f) findet man leicht, daß 
Y ex ex 
E(t +) Eg(t) — (ft — w) 


Dann bleibt mit (33) 


N 1 4 ne? 1 
evr" 4 VN 2Mo Wiw) "\2 r Q W(w) 


1 | 


lee, m) G(x) [ Votr.) (0) 1 2, w* 
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Zur Berechnung der verbleibenden Terme bilden wir z. B D: 
<f 

1 
B 
k 
st 
a 
d 
| 
S 
Ä 
(50) k 
d 

u 
is 
t I 

ex | 
| | 
4 

ars 

| 
5 
Bs (51) | 
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Damit gehen wir in (31) ein und finden mit (50) 


1 b 52 
w 
die Diskussion des Resultats = N 
ent- Der Vergleich des Matrixelements (52) mit dem entsprechenden Resultat 
Bardeens zeigt zunächst, daß es hinsichtlich der Übergangswahrscheinlich- 
keiten keinen Unterschied macht, ob man in nullter Näherung Bloch -Zu- 
(49) stände benutzt und den Wechselwirkungsoperator linear in $ ausdrückt, oder 


aber von Wellenfunktionen ausgeht, welche mit den Ionen ‚‚mitgehen‘‘, derart, 
daß sie in der Umgebung eines Ions im wesentlichen ihre alten Werte behalten, 
en: und dabei die Wechselwirkung linear in den Deformationen &, „ ausdrückt. Diese 
Situation wird verständlich vor dem Hintergrund der adiabatischen Näherung, 
welche schon dem Ansatz (6) zugrunde liegt. Wenn die Elektronen den Ionen- 
lagen adiabatisch folgen, so kann offenbar die Verschiebung der Ionen ohne 
Änderung ihrer relativen Lage zueinander keine Übergänge bewirken, diese 
50) können erst durch Veränderung der relativen Konfiguration, also Verzerrung 
des Gitters, verursacht werden. Die Verzerrung des Gitters zieht eine Defor- 
mation der Wigner-Seitz-Zelle nach sich, welche Anlaß zu Übergängen 
gibt, während durch eine Translation der Zelle ohne Deformation nur eine 
adiabatische Veränderung der Wellenfunktion ohne Übergänge erfolgt. Dar- 
aus ist zu ersehen, daß für die Übergänge nicht 8, sondern Eoo maßgebend sein 
muß. Da sich die Funktionen yr und Pr = e”iS yr nur um einen Term O (e,,) 
unterscheiden, weil yr das Produkt aus zwei periodischen Funktionen ist, ist es 
bei einer Rechnung, in der die Störung im Hamiltonoperator linear in &,. 
ist, gleichgültig, ob man mit yr oder Y; die Matrixelemente berechnet, wenn 
man diese nur linear in Eo haben will. 

Was die Ankopplung der transversalen Gitterschwingungen tiber Normal- 
prozesse betrifft, so ist sie, bei Beachtung des Energiesatzes, von der Form 


h2 
+ w {Vj + trans. ~ h o-(ep[1 (te) [1 


m* m* 


und verschwindet für m = m*. Da die atte! Ye der longitudinalen Phononen 
9 
von der Größenordnung 3 w Ep ist und * „ (m f) (fe) » 2Er-w, wird sie 


m 


vergleichbar mit der der transversalen, wenn 1— = > 0,3 ist. Mit dem 


m r r 
Wert 1,29 von welcher etwa bei Hunter und Nabarro®) für Na!) ange- 
m 


U 

geben ist, kommt man also bereits an den Fall heran, in dem longitudinale 
und transversale Phononen konkurrieren könnten. Dabei hat man jedoch zu 
beachten, daß nur diejenigen Elektronen mit den transversalen Gitterschwin- 
1) gungen in Wechselwirkung stehen, für die sowohl (tw, f) + 0 als auch (e,D + 0 


12) Dieser Wert ist sicher zu groß, so daß wahrscheinlich die transversale Kopplung 
für Normalprozesse praktisch verschwindet in diesem Falle. 


365 
= 
4 
= 
3 


ist, während die übliche Kopplung an die longitudinalen Phononen unabhängig denn 
von f ist. Dadurch wird der Anteil der transversalen Phononen an Normal. 
prozessen stark reduziert. Die Frage nach dem Einfluß dieser Wechselwirkung 
auf die Transporteffekt ein Metallen soll späteren Untersuchungen vorbehalten 
bleiben, wobei auch die Bedeutung der hier ganz außer acht gelassenen Um- 


nd 
klapp-Prozesse betrachtet werden muß. ’ 
AnschlieBend sei Herrn Dipl. Phys. H. Puff fiir das in zahlreichen Diskus- Reo 
sionen bekundete Interesse und fiir wertvolle Hinweise herzlich gedankt. a 
nd 
Anhang 
Zum Beweis der Formel (41) merken wir zunächst an, daß wegen (37) und 
der Hermitizität des Operators bz, Oz, 
f) = De. (E) = Dgalk), 
Deck) = Dy o(\€!) Do o(k). Wis: 


Nun mögen durch die zyklischen Indizes x, ß, y die drei orthogonalen Koor- 
dinatenachsen bezeichnet sein und es sei S‘ die Spiegelung an der a B- Ebene, 


also 


— 
0 0 -1 
Es ist [8] = 8%) und || S®|| = — 1. Da die Spiegelungen S® alle zur 


Punktgruppe der kubischen Gitter gehören, gilt 


us (Sr t) = ut(t). 


Mit 


haben wir daher für 9 + a, (0, 0, t zyklisch) 
82 ur) (t’) 


(r) 
Dealt = u*(r) y(t’) - x’, || dx 


3 Dg o(k) Dee (k)- 
Nun sei 7” eine Drehung um die t-Achse um den Winkel >. Es ist || 7 || = 
+1 und mit "= 


02 (t’) 
Ox, Ox, 


Deo (t) = 2 The: || 

17} 


Perg 
ud 


02 Up(r) y(t’) 


0x’? 
0 


I 


| u 


d.h 4 
| 
= 
+ 


ınd 
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denn es ist z. B. 


= 
und allgemein = 


(Tre = (0,0,r zyklisch) 
Doo lk) = Dso(|TE|) = Die k) 
und daher 


sf 02 1 P 


2, 
= (ur Aug + up Auf) dr = — 


Berlin, Physikalisch-technisches Institut der Deutschen Akademie der 
Wissenschaften. 


Bei der Redaktion eingegangen an 6. Juli 1960. 
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rperiodischen Potential 
Vo(t) = Vo(t + 9) (1) 


stehendes Elektron und fragen nach Fin h eine Gitterdeformation ausgelösten 
Übergängen zwischen seinen stationären (Bloch-) Zuständen. 
Das ungestörte Problem wird beschrieben durch den Hamilton-Operator 


H, = Ha + Hon 


Wir betrachten ein unter dem EinfluB eines gitte 


1 

Ha p? + J g(t), Hn > op, ; (bw, ; bi, + ; dw, ;) (2) 

® Dabei ist wegen der Separierbarkeit des Problems 


Die mit 

A 

= Ey. pe (t) = ug (et), (t) = (t+ 9), = 
i i 1 i 


w, 7 


(4) 


Wird die Gitterdeformation durch H,, beschrieben, so ist die sekundliche 
Übergangswahrscheinlichkeit in erster Näherung gegeben durch 


(Gj) \ N | (/) \ 
m,3 (5) 
ÖölEr — Er + S — Ny, j)h wy, ; 
(Da uns nur Übergänge zwischen Elektronenzuständen interessieren, summieren 


wir gleich über die Endzustände des Phononenfeldes). Machen wir nun für 


das Potential des deformierten Gitters den bekannten — 


h 1/2 » ew.j i + i 
va) > zei + bi, zer), (6) 
iD, 


weite 


Mitte 


so di 


er 
Von Heinz Puff 
ru 
Im f 
| 
gilt. 
| = 
Berii 
Weis 
mit 
ni 
Yr 
Ce 


entials 
(1) 
lösten 


erator 


(2) 


dliche 


nieren 
in für 


(6) 
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/ /2 
a, (u(r), grad Vo) -6 V (7) 


ei (t,t) + bw. e~ (Cw, grad 


Setzt man diesen Ausdruck in (5) ein und führt die Matrixelemente bezüglich 
der Phononen-Operatoren aus, so ergibt sich in bekannter Weise: I 
1 A~ Ge 


n 


wrt (tw) = 
(em,;, grad Vo)| ölEr — Er— (8) 
+ (rw; + 1) (ew, Vo)|yr > PölEr — Et + hor, 
Im folgenden benutzen wir, daß für jede (skalare) Funktion G des Ortsoperators 
t die Relation 


[p, @(r)] = grad, G(r) 


gilt. Damit ergibt sich: 
i + >> 
= (pr je* i(tw, Cw, „grad = + le i(w,r) P). V. ol! ve 
{pr le‘ p)| Voyr>— <Voyr le‘ i(w,t) (ew, js p)! Yr 
weiter wegen wr = Et yr: 
Mittels (9) folgt lei 
littels (9) folgt leicht 


[p?, et = + 2 Alm, p) et — 


so daß wir schließlich erhalten: oe ye ur 
My: = — Er — je‘ itw,r (Ew, p)| yp } 


ens der expliziten Form der Bloch -Zustände liefert in gleicher 
else: 
Yr le (eu; p)| yr> (ur jet (ew, 5, P)| ur) 

+ h(ew,j; f) Cup je‘ wet uy » 

(pr p) (ery, ye> = Cur (rv. p) (Cw. 
+ A(w, + w (ur fe (ey 5, p)| 
+h (ew, Sur (m, p)| ur) 
+ RR? (Cw, „Di, +8 cur | | 
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a 
¢ 4 
5 | 
(3) 
(4) 
| 
(5) 


und damit insgesamt nach Reduktion aller Matrixelemente auf Integral 
über eine Elementarzelle: 


h? mw? ‘ 
2m 


5 
{[# — By + 


+ = (Cw, ;, ia ut (r) (tw, grad) ur (t) 


FAR ub (m, grad) grad) m). 


Ay 


(ur = grad, By — ht 


jet (ey ;, grad Vo)| = + 


f) [ (ww, grad; Er) — | 
“39 


* (Cw, gradı Er) — = dr uf (x) (tv, grad) (fi, ;, grad) ur 


und mittels (8) ergibt sich schließlich: 


1 


i | he 
wre lo-0 = WG 2 | (ew,;j- [ (w, grad; Er) — — (m, 


m 
bi [= (tv, f) — how, i| (Cw, j, grad; E 


m 


+ ; + 1) 6(By — Er + Ör-t,-w}- 


Sarutte) (v0, grad) grad)ur(e) | (Er 


h? 9 | 
% uF (t) (ew,; grad) ur(t) 


(11) 


in quadratischer Terme sowie mit der Näherung 
er & ur geht (11) unter Berücksichtigung der Relation iii 
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Die räumliche Intensitätsverteilung der Keilinterferenzen 
Von Günter Schulz und Günter Minkwitz 

Mit 7 Abbildungen 


Inhaltsübersicht 


Die räumliche Struktur der Keilinterferenzen, wie sie z. B. in Michelson- 
Interferometern auftreten, wird berechnet und dargestellt. Die Ergebnisse 
liefern die Intensität als Ortsfunktion im ganzen Raum. Sie erscheint im 
wesentlichen durch den Kontrast und eine die Streifenlage bestimmende 
Größe eindeutig charakterisiert. Dabei werden auch spezielle Verknüpfungen 
zwischen den Keilinterferenzen und den Beugungserscheinungen bei aber- 
rationsfreier Abbildung aufgezeigt und ausgewertet. 


1. Einleitung 


Keilinterferenzen treten in Michelson- und Mach-Zehnder-Interfero- 
metern, in Interferenzmikroskopen und anderen optischen Anordnungen auf. 
Ihre Behandlung bildet die Grundlage fiir die Erfassung der Interferenzen 
gleicher Dicke oder Fizeau -Streifen. Man hat immer wieder versucht, sie 
theoretisch möglichst vollständig zu erkennen. Aus der zahlreichen Literatur 
seien nur einige Beispiele aus der letzten Zeit genannt. 

Schulz!) gibt eine lineare Näherung für die räumliche Struktur der Keil- 
interferenzen, verbunden mit einer exakteren Bestimmung der Intensität in 
der Ebene ihrer größten Schärfe. Hopkins?) wendet Methoden der Kohä- 
renztheorie an und untersucht damit Erscheinungen in Zweistrahlinterfero- 
metern. Hansen und Kinder?) diskutieren ebenfalls den Kontrast der 
Fizeau -Streifen. Landwehr) knüpft an frühere Untersuchungen von Igna- 
towsky®) an und gibt eine Darstellung von Intensität, Sichtbarkeit (Kontrast) 
und Streifenlage in der Ebene der einen Keilfläche. Lohmann’) gibt eine 
quadratische Näherung für den räumlichen Kontrastverlauf der Keilinter- 
ferenzen. 

Im folgenden soll nun eine umfassendere Darstellung der räumlichen 
Struktur der Keilinterferenzen gegeben werden. In der Entwicklung ihrer 


1) G. Schulz, Ann. Physik (6) 14, 177 (1954). A; 

*) H.H. Hopkins, J. Opt. Soc. Amer. 47, 508 (1957). 

3) G. Hansen u. W. Kinder, Optik 15, 560 (1958). 

‘) R. Landwehr, Z. Naturforschg. 18a, 58 (1958). 

5) R. Landwehr, Opt. Acta 6, 52 (1959). Vgl. auch das vorhergehende Zitat*). 

*) W. v. Ignatow sky, Einfluß der Form und der Lage der Lichtquelle bei den 
Messungen mit dem Interferenz-Komparator nach Kösters, Moskau 1935. 

7) E. Ingelstam, Problems related to the accurate interpretation of microinterfero- 
grams, Symposium on Interferometry Teddington 1959, paper 3—2, p 4. 


$ 
| 
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Form spielen bei kreisförmiger Blende die Lommelschen 
| Funktionen ®)?) eine Rolle, wie insbesondere auch bei Landwehr?) ersichtlich 
ist. Im Gegensatz zu den Arbeiten von Ignatowsky und Landwehr jedoch 
steht bei der vorliegenden Abhandlung die gesamte räumliche Intensitatsver. 
teilung im Mittelpunkt der Darstellung. Das erfordert nun zweckmäßig eine 
andere Methodik und dabei insbesondere von vornherein ein anderes Koordi- 
natensystem, nämlich dasjenige, dessen Koordinatenebenen im Raum Sym- 
metrieebenen der Interferenzerscheinung sind. Wir schließen uns daher hier 
an die Arbeit!) an. 


2. Optische Gegebenheiten 


Wir gehen von folgenden Gegebenheiten aus. Durch eine optische Anord- 

nung werden im Unendlichen zwei miteinander kohärente Bilder einer Licht- 

quelle erzeugt, die durch eine reine Drehung im Endlichen zur Deckung 
gebracht werden könnten. Als Beispiel einer solchen ar 

i<”s optischen Anordnung ist in Abb. 1 ein virtueller Keil in we on 3 

einer Michelson-Schaltunggewählt. DortwirddieLicht- 

ip oso Q in die Ebene der Aperturblende B abgeb’ldet. 

_ B liegt in der hinteren Brennebene der Linse Z, durch 

A das Licht über den teildurchlässigen Spiegel 7 auf 

die beiden Spiegel R und $ fällt. Diese reflektieren das 

Licht über 7 nach oben, wo es zur Beobachtung gelangt. 

i 8’ ist das durch T entworfene Bild von S. S’ bildet mit 

Reinen virtuellen Luftkeil mit der Keilkante ce und dem 

 Keilwinkel 9. Der wirksame Aperturkegel!)!%) sei ein 


Lichte. 


Beobacht 


Abb. 2 
Optischer Raum des 
Keils; Keilkante c 
senkrecht auf der 
Zeichenebene (Keil- 
flächen selbst sind hier 

-1%7 virt Kei nicht eingezeichnet, da 
... ud sie für die Lage der In- 
J terferenzerscheinung 

£ c nicht unbedingt cha- 
Abb. 1. Virtueller Keil in einer Michelson-Schaltung rakteristisch sind) 


R 


Kreiskegel; die wirksame Apertur u nehmen wir hier klein gegen 1 an. 
Aberrationen, Beugung und besondere Polarisationseffekte werden vernach- 
lässigt. Das Licht wird als monochromatisch angenommen (Wellenlänge 3). 

Die dabei entstehenden Keilinterferenzen betrachten wir im optischen 
Raum jenes Keils, siehe Abb. 2. Unser Koordinatensystem orientieren wir dabei 


nach der Keilkante ce und der senkrecht dazu angenommenen mittleren Licht- 
Boe %) E. Lommel, Abh. d. math.-phys. Cl. d. kgl. Bayr. Akad. d. Wiss. 15, 229 (1886). 
ie m) 9) G.N. Watson, A Treatise on the — of Bessel Functions, Second Edition 
Cambridge 1958, S. 537. 
bat 7 10, G, Schulz, Optik 16, 280 (1959). 
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richtung!)). Diese halbiert in jedem Punkt P (p,h) den Winkel zwischen 
denjenigen kohärenten Strahlen k’ und k”, die durch den Mittelpunkt der 
Aperturblende gegangen sind. c selbst ist die Koordinatenachse p = 0, h = 0. 
Die dritte Koordinate, parallel zu c, wird nicht benötigt, denn die Interferenz- 
erscheinung erweist sich innerhalb unserer Voraussetzungen als von ihr un- 
abhängig. Wir nehmen dabei an, daß für jeden betrachteten Punkt P alle 
Strahlen innerhalb der beiden schraffierten Kegel (die sich auch beliebig 
durchdringen dürfen) noch zur Beobachtung gelangen. Die Lichtamplitude 
soll unabhängig von P und von der Strahlrichtung sein, wobei nur ein kon- 
stantes Amplitudenverhältnis zwischen den kohärenten Strahlen zugelassen ist. 


3. Ableitung der Formeln 


Unter den obigen Voraussetzungen erhält man für die Intensität /, nach 
entsprechender Normierung, folgenden Integralausdruck!): 


u 22 


- f J cos (s cos p +¢tsing cos x)singdpdz. (1) 


I — COS u) 
Hier bedeuten: A > 1 eine ortsunabhängige Konstante, die z. B. vom Ver- 
spiegelungsgrad von T in Abb. 1 abhängt und optimal 1 ist; s = 22 p/C und 
t= 2ah/C Ortskoordinaten, die zu den Ortskoordinaten p bzw. h propor- 
tional sind; dabei ist . ‘enue der Streifenabstand, bis auf Korrektur- 


faktoren, gleich C = — 10), also im s-t-MaBstab 2 2. Nach Anwendung 


sin d 
des PERF OEEER.; auf den Ausdruck cos (s cos g + ¢ sin @ cos x) und 
Ersetzung von cos (sing cos x) durch cos (¢ sin p sin x) für die Integration 


über eine volle x-Periode erhält man aus (1): u np. = 


I=4+ cos (s cos cos (f sin @ sin sin g dq dx, 


I=4+ J cos (9.605 ) Jo (t sin sing dp 


1 — cos u 


mit der Besselfunktion nullter Ordnung AEN 


J,(y) = ee | cos wa sin in x) dx. 


u 


9 
f=A+— J cos [s(1 — g/2)| Joltp) p de. 
z COS 8 J cos (s 92/2) p) dg | 
74 + sins | sin(s g?/2) Jg(tq) 
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Die letzten beiden Integrale sind bereits von Lommel’) bestimmt worden, 
so daß wir schreiben können: 


: I=A + e0ss(U, cos — + U, sin | 
(4) 
sin s(U, sin > 2 COS 2) | 
ey mit den Lommelschen Funktionen?) U, und U,: 
Uy = 0, (w,2)= 
ist dabei die Besselfunktion n-ter Ordnung, während 


tat it ist, vgl. auch den Text hinter (1). Aus (4) erhält man | weiter: 


I=4+ U, cos (s = 5) U, sin (+ (10) 
Führt man hier 


is ein, dann kann man in (10) das Additionstheorem anwenden und erhält da- 
durch 


= I=4A+Kcos(s—®) (14) 


Speziell für w — 0 wird nach (5) bzw. (6): 


Ents 


hier! 


eilt 
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(11) 
(12) 
(13) 
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(14) 


(15) 


@. Schulz u. G. Minkwitz: Intensitätsverteilung der Keilinterferenzen 375 


also darf man En für (10) bzw. (13) schreiben: 


= ı(2) (w=s=p=0), (16) 


: 
= Ji (2) (w=s=p=0).— (mM 
Entsprechend wird für z — 0: bl 
hiermit folgt aus (10) nach des 
4 12 
I = A+—~- cos (s — ) (2 =t=h=0), (18) 


während sich aus (13) nach entsprechender Umformung der Winkelfunktionen 
ergibt: = 


(19) 


Definiert man nun 


| (20) 
so kann man die Struktur der Linien m = const in der w-z-, s-t- bzw. p-h- 
Ebene aus ® bestimmen. Dies ist im allgemeinen, wegen der meist schwachen 
s-Abhängigkeit von ®, recht einfach (siehe Abschn. 4). Durchführbar ist es 
jedoch auch bei stärkerer s-Abhängigkeit, die in der Umgebung singulärer 
Punkte auftritt. Als Beispiel bringen wir kurz die Bestimmung der Neigung 
der Linie m = 1 für z = a, (erste Nullstelle von J, (z)); der zugehörige w-Wert 
sei w,. Da man zeigen kann, daß für z = a, 


| 
gilt, ergibt sich für w,: 
er 
wofür wir wegen der vorausgesetzten Kleinheit von u den Wert ? = 7 u? setzen. 


Die gesuchte Neigung im Punkte (w,, %) —_ dann aus: 


1) Die gleiche Form bereits bei Hopkins?), Gleichung (42). 
12) Diese Gleichung und ihre Darstellung bereits in !), und zwar exakt auch für größere 
u; siehe dort Gleichung (21) und Abb. 4. 


1) K entspricht in 1), 


(5) 
(6) 
(7) 
K 1 
= 
- en 
7 4 
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wobei hier m als Funktion von w und z mit dem Parameter u? aufzufassen ist. 
Geht man zur Grenze u? — 0 über, erhält man: 


3 

az a 

dw a: 2a, 1,230 (m 1, Oy 3,83) (21) 


= —1,230u(m@ =1,1=% = ). (22) 


dp ds 20, 


4. Ergebnisse für die Keilinterferenzen 


Aus (14) entnehmen wir, daß die Intensität / als Ortsfunktion im wesent- 
_ lichen aus einer modulierten eos-Funktion besteht. Die cos-Funktion selbst 
ist für das Vorhandensein - Interferenzstreifen verantwortlich: Überall 
dort, wo nach (20) m = 0.1, 2,... ist, liegt ein Intensitätsmaximum; d.h. 
dort geht der en der Ordnung m hindurch. m ist dabei eine 
„mittlere Ordnung‘, die durch Mittelbildung über die Wirkungen der einzelnen 
inkohärenten Elemente der Aperturblende entstanden ist. Die genaue Lage 
der Streifen wird durch die Korrekturgröße ® bestimmt; dabei kann man, 
nach (20), ® = s— 2nm als Abweichung eines Streifens gegenüber seiner 
unkorrigierten Lage s = 22 m auffassen. Der Kontrast oder die Sichtbarkeit 
= beträgt: 

max min 

X 


~ (A+ A’ (23) 


; 
Abb. 3. K als von w und z 
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wobei A > 1 die hinter (1) erläuterte Konstante ist und hier nicht näher 
interessiert. K und ® hängen nach (13) und (15) vermittels (5) bis (7) und (11) 
bis (12) nur von w und z ab. Insgesamt bedeutet das nach (8) und (9), daß die 
gesamte Interferenzerscheinung — für alle Wellenlängen 7, Keilwinkel 9. 
kleinen wirksamen Aperturen « und im ganzen Raum (p,h) — quantitativ 
vollständig und übersichtlich bestimmt ist, wenn man K und ® in Ahhängig- 


keit von w und z kennt. Diese Abhängigkeiten sind nun nach den erwähnten 


= = 
4 = 


2 4 6 6 10 12 


Abb. 4. ® als Funktion von w und z; zur besseren Darstellbarkeit ist — ® aufgetragen; 
der — ®-Wert der gezeichneten Grundebene beträgt — 7 


Gleichungen in den Abbildungen 3 bis 5 veranschaulicht. w und z kann man 
dabei als allgemeine Ortskoordinaten auffassen, die nach (8) und (9%) zu den 
ursprünglichen Ortskoordinaten s und ¢ bzw. p und h verschieden proportional 
sind, je nach der Größe von u. 

Einige Spezialfälle der genannten Abhängigkeiten sind folgende. w= 0 
bedeutet die Ebene durch die Keilkante parallel zur mittleren Lichtrichtung. 
Dort sind Intensitäts- bzw. Kontrastverlauf parallel zur mittleren Lichtrich- 
tung nach (16) bzw. (17) im wesentlichen durch die Besselfunktion 1. Ord- 
nung bestimmt. So erhält man für A den in Abb. 3 über der z-Achse aufge- 
tragenen Verlauf, wobei \,. 9, . - . die Nullstellen jener Besselfunktion sind. 
>= 0 bedeutet die Ebene durch die Keilkante senkrecht zur mittleren 
Liehtriehtung, d. h. die Schärfeebene (hier zeigen die Interferenzen im wesent- 
lichen die größte Schärfe). In dieser Ebene sind Intensitäts- bzw. Kontrastver- 
lauf senkrecht zur Keilkante aus (18) bzw. (19) erhältlich. Für A führt dies 
zu dem in Abb. 3 über der w-Achse in perspektivischer Sicht aufgetragenen 
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Verlauf, welcher unverzerrt die in Abb. 6 gezeichnete Form hat. — Für kleine 
w und z erhält man als Linien mit gleichem K in der w-z-Ebene Ellipsen, 
deren Achsenverhältnis mit dem in einer von Lohmann‘) angegebenen Nähe. 


mittt.| Lichtr. 


Schnitt 
mit der 
Schörfe - 
ebene 


Abb. 5. K in der w-z-Ebene, d.h. in einer beliebigen Ebene senkrecht zur Keilkante ¢; 
die K-Werte sind bei den verschiedenen Schattierungen, in Bildmitte beginnend: 
1 = K > 0,866 (weiß), 0,316 > K > 0,224, 
0,866 > K > 0,707, 0,224 > K > 0,158, 
0,707 K > 0,500, 0,158 > K > 0,100, 
0,500 > K > 0,316, 0,100 > K>0 (schwarz) 
rung übereinstimmt, wie sich zeigen läßt. — Ferner ergibt eine andere Speziali- 


sierung unserer Ergebnisse einen weiteren bereits gerechneten Fall: Bei Land- 
wehr?) sind, in unserer Bezeichnungsweise ausgedrückt, X und ® auf einer 
schwach geneigten Keilfläche <1 in Abhängigkeit von einer apertur- 
bedingten Größe ange- 
geben; die Umrechnung 
. . %Y . 
zeigt auch hier Überein- 
stimmung, nachdem eine 
weitgehende formale Ana- 
logie dies bereits vorher 
vermuten läßt. 
Aus ® läßt sich die 


genaue Lage der Inter- 


0 | | 
0 47 67 eI ferenzstreifen bestimmen. 
Abb. 6. K für z = 0 als Funktion von w allgemeiner, die genaue 

Lage der Flächen konstan- 

ter mittlerer Ordnung m im Raum. — Z. B. ist nach Abb. 4 in der Schärfe- 


ebene z = 0 ® abschnittsweise eine lineare Funktion von w, und jeweils in 
einer gewissen Umgebung davon ist diese Linearität wenigstens noch ange- 


nähert erfüllt. Im ersten Abschnitt 0 <w <4r ist dabei ® 1 87: 
so daß wir (20) in der Form schreiben können: zu 
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Der Streifenabstand im s-t-MaBstab beträgt also hier nicht genau 2 2, sondern 
erscheint um den bereits von Ignatowsky®) berechneten Korrekturfaktor | u 
|+ u?/4 vergrößert. — Wo ® dagegen nicht mehr als lineare Funktion von 
rund z und damit von s und ¢ gelten darf, läßt sich ein derartiger Korrektur- 
faktor nicht mehr angeben; die Streifen haben örtlich etwas verschiedenen 
Abstand voneinander. Die Lage jeder Fläche 

nm = const im Raum ist dabei immer noch durch i 

Gleichung (20) gegeben, welche fiir jedes feste 7 verschw. 
i eine Gleichung zwischen s und ¢ mit dem Para- FF Kontrast 
meter u darstellt. Ihre Lösung wird dadurch er-- 
kichtert, daß ® an den meisten Stellen zwar 
och wesentlich von w und z, aber nur noch schwach 
von s und £ abhängt. Man braucht dann in der 


mit \lichtr. 
Gleichung 
s=-2rzm+P® 
® selbst nicht genau an der Stelle s zu nehmen, verschw 


sondern darf den ®-Wert der Nachbarstelle 22 m 
einsetzen, wenn ® nicht zu groß ist; so erhält 
man s explizit als schwach t-abhängige Funktion. | 

- Ist auch dieses Näherungsverfahren nicht mehr Abb. 7. Schnitt | 
- Abb. 7. Schnitt | ¢ durch 
möglich, dann stehen andere zur Verfügung. Als die Flächen ganzzahliger 
Beispiel verweisen wir auf das zu (22) führende Ver- mittlerer Ordnung min der 
fahren, welches, zusammen mit dem zuletzt ge- Umgebung der Stellen ver- 
schilderten, zu dem in Abb. 7 dargestellten Bild ne eae Kontrastes 
ei Mm 0; das Bild gilt 
führt. Man erkennt (vorausgesetzt, daß derartig jm p-A- und im s- Malistab 
kleine Kontraste noch wahrnehmbar sind): Fokus- für wx 

siert man, von der Keilkante e ausgehend, parallel 

zur mittleren Lichtrichtung, so verschwindet zu einem bestimmten Zeitpunkt 
der 0. Interferenzstreifen, zusammen mit dem Kontrast, exakt. wobei die — 
beiden Nachbarstreifen ziemlich schnell auf den normalen Streifenabstand. 
zusammenrücken. 


er Kontrast 


Verknüpfungen mit Beugungserscheinungen 
Bei der oben geschilderten Struktur der Keilinterferenzen hatten wir 
Beugungseffekte vernachlässigt. Interessant ist nun, daß sich zwischen dieser _ 
Struktur ohne Beugung und gewissen anderen, gerade durch Beugung beding- 
. . . Y . . . 
ten Erscheinungen sehr weitreichende formale Übereinstimmungen aufzeigen _ 
lassen. 
. . . 
Daß zwischen beiden Gebieten Entsprechungen bestehen, wird z. B. auch 
dadurch nahegelegt, daß die in (3) stehenden Integrale auch bei den von. 
Lommel) behandelten Beugungserscheinungen auftreten. Diese Entspre- 
chungen gehen so weit, daß z. B. die in-(11) bis (13) auftretenden Ausdrücke 
) » 
U, und = U, bereits in der Lommelschen Arbeit tabelliert sind, ebenso 
das Quadrat von K. So wurde fiir die Abbildungen 3 und 4 fast keine zusätz- 
liche numerische Arbeit nötig. 


Zur Darstellung von Abb. 5 dagegen reichen die Lommelschen Rech- 
nungen nicht aus. Dafür kann jedoch eine Verknüpfung mit einem von | 


| 
kleine 
ipsen, | 
Nähe. | 
| 
1 
> 
: 
¢ 


380  Annalen der Physik. 7. Folge. Band 7. 1961 


Berek") behandelten und auch von Picht!?) gebrachten Beugungsfall auf. 
gezeigt werden. Es handelt sich dort um die aberrationsfreie Abbildung eine 
Punktes unter Berücksichtigung der Beugung. Läßt man nun folgende Ent. 
sprechungen zwischen unseren Größen und denen bei Picht gelten: 


— Wellenlänge, u 
— Abstand von der Bildebene. a 
— Abstand von der optischen Achse, 


— Apertur, 

l 2, 

dann folgt, daraus wie ein Vergleich der Formeln zeigt, die Entsprechung 

Wellenlänge 
K — Intensität. 
a - Apertur? 
Man kann also die Bereksche Abbildung über die Intensitätsverteilung in 
. der Umgebung des Bildpunktes direkt für unsere K-Verteilung übernehmen, 
u. wenn man die entsprechenden Größenumdeutungen und Maßstabsänderungen 
@ . . - . . 

: vornimmt. Dies ist in Abb. 5 geschehen. — Bemerkenswert ist hier noch, daß 


a zwar die Entsprechungen bei A und « im „richtigen‘‘ Sinne vorhanden sind, 


: daß jedoch die Koordinatenrichtungen gerade vertauscht erscheinen, d.h. 
. . . . . 
der parallel zu A verlaufenden mittleren Lichtrichtung entspricht nicht die 


optische Achse, sondern gerade die Richtung senkrecht dazu. 


14) M. Berek, Z. Physik 40, 420 (1927). ee 

15) J. Picht, Optische Abbildung, Braunschweig 1931, § 21. 


Berlin-Adlershof, Institut für Optik und Spektroskopie der Deutschen 


Akademie der Wissenschaften zu Berlin. 
Bei der Redaktion eingegangen am 27. Juli 1960. ny Oe 
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Zur Realstruktur von CdS-Einkristallen 
Von W. Méhling und H. Thiel 


Mit 6 Abbildungen 7 


Inhaltsübersicht 


Die Winkeldesorientierungen der (0110)- und der (0001)-Ebenen von 30 
aus der Dampfphase gezüchteten CdS-Blättchen wurden bestimmt. Ein Zu- 
ammenhang zwischen Zahl und Stärke der Rippen und der Desorientierung 
wurde nicht gefunden. Verschiedene Kristalle waren gekrümmt. Es wird . 
angenommen, daß diese Krümmungen durch Bildung überse hüssiger Stufen- 
versetzungen eines Vorzeichens während des W achstums bzw. während der 
ersten Phase der Abkühlung entstehen. Das Aufkleben der Blättchenkristalle 
verspannt das Gitter im allgemeinen beträchtlich. Ein eventueller Einfluß 
auf die Leitfähigkeit, selbst bei Beanspruchung bis zur Bruchgrenze, liegt 
unterhalb 19%. 


D 


Viele physikalische Eigenschaften werden durch Kristallbaufehler ent- 
scheidend beeinflußt!)?). Linienförmige Fehlstellen, sogenannte Versetzungen, 
sind vor allem für die plastischen Eigenschaften der Kristalle verantwortlich 
und wurden bisher vorwiegend am Silizium und Germanium untersucht !)?). 
Mit analogen Untersuchungen ist neuerdings auch an C'dS-Einkristallen be- 
sonnen worden 


Meßanordnung 


Für unsere Untersuchungen wählten wir die Methode nach Guinier und 
Tennevin®), im Prinzip ein fokussierendes Laue-Verfahren (Abb. 1), mit der 
man zwar nicht einzelne Versetzungen nachweisen, wohl aber mit großer 
Genauigkeit die Paralleliti ät der Netzebene n bili ätte :henförmiger Kristalle kon- 
trollieren kann. 


1) E. Billig, Growth and Defects of Se mirseelnsten Crystals. Halbleiter und Phos- 
phore. gsbe richt Garmisch 1956. 

2), P. Haasen u. A. Seeger, Plastische Verformung von Halbleitern. Vortrag Phy- 
sikertagung Heidelberg 1957. 

4) D.C. Reynolds u. L.C. Greene, Crystal Growth Mechanism in CdS. J. appl. 
Phys. 29, 559 (1958). 

4) E. Votava, Versetzungen in CdS-Kristallen. Z. Naturforschg. 13a, 542 (1958). 

5) A. Guinieru. J. Tennevin, Sur Deux Variantes de la Méthode de Laue et leurs 
Applications. Acta Cryst. 2, 133 (1949). 
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Das von einem strichförmigen Brennfleck S einer Röntgenröhre mit Wolf. 
ram-Antikathode ausgehende schwach divergente Strahlenbiischel wird in 
Kristall A abgebeugt und in F fokussiert, da die virtuelle Lichtquelle S’, de 
Kristall A und der Reflex F auf dem Fokussierungskreis liegen. Bildet de 
Mittelstrah! SO des Strahlenbüschels mit der Kristallnormalen den Winkel ı 
dann ist nach der Reflexion in O de 
Winkel zwischen Primärstrahl und 
reflektiertem Strahl 2x. Bei einen 
Öffnungswinkel 22 des Strahlen. 
büschels wird in A der Auftreffwinkel 
auf die Netzebenen um e kleiner und 
in Bum e größer sein als inO. Da es 
sich um weiße Strahlung handelt, wird 
bei nicht zu großem e auch hier Strah. 
lung reflektiert, jedoch unter einen 
entsprechend kleineren bzw. größerer 
Winkel. Das vor der Reflexion diver. 
gente Strahlenbüschel wird nachher 
also konvergieren®). Längsausdehnung 
des Brennflecks und der Netzebenen 
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Abb. 1. Laue-Verfah- 


Das fokussierende 
ren nach Guinier und Tennevin (sche- 


rohre, senkrecht auf der Zeichenebene; 
K = CdS-Kristall;a = Richtung der a-Ach- 
se in der Zeichenebene (c steht senkrecht auf 
a der Zeichenebene); F fokussierter Re- 
flex; S’ = der zu S in bezug auf K sym- 


q 
1 
; matisch). S= Strichbrennfleck der Röntgen- 


metrische Punkt auf dem Fokussierungs- 


stehen senkrecht auf der Zeichenebene. 


kreis; M — Mittelpunkt des ee aig Die Breite des Reflexes im Fokussie- 
kreises; 2 ¢ = Offnungswinkel des Rönt rungspunkt F ist bis auf gewisse 
ö genstrahlenbüschels; x = Winkel zwischen ma leich d B fleck 
Mittelstrahl des Büschels und den beu- orrekturen gleich der brenntleck- 


genden Netzebenen 


breite. Sind die Netzebenen entgegen 


der ursprünglichen Annahme nicht 

parallel, so ergibt sich eine Reflexverbreiterung, aus der man auf die Desorien- 
tierungen zwischen den Kristallbereichen schließen kann. 

Das Auflösungsvermögen wächst mit der Fokussierungsentfernung OF 

Bei dem verwendeten Aufbau betrugen die Abstände SK und KF jeweil 

100 em. Winkeldesorientierungen bis zu + 15 see konnten bestimmt werden. 


3. MeBergebnisse 

Die von uns untersuchten CdS-Einkristalle sind blättchenförmig (8— 10 mm 
lang, 2—5 mm breit und 0,05—0,1 mm dick) aus der Dampfphase gezüchte 
und haben Wurtzitstruktur’); kubische Anteile oder Lagenfehlordnung 
konnten bisher nicht beobachtet werden®). Die Kristallblättchen zeigen ein 
mehr oder weniger große Anzahl von Streifen oder Rippen, die parallel zuein- 
ander angeordnet sind und senkrecht oder schräg auf der Längskante de 
Kristalls stehen. Diese Riffelung verläuft parallel zur c-Achse, während di 
a-Achse mit der Blattchennormalen zusammenfällt®). Untersucht wurden je 
weils die Reflexe der (0110)- und die der (0001)-Ebenen. Die Reflexbreite! 

®) Der nicht abgebeugte Teil der Strahlung bleibt selbstverständlich divergent. 

7) F,Ullrich u. W. Zachariassen, Über die Struktur des a- und f-CdS. 
Z. Kristallogr. 62, 260 (1925). 

8) W.Schweinberger, Röntgenographische Untersuchungen an (dS-Kristalleı 
Diplomarbeit Berlin 1955. 
9) H. Gobrecht u. A. 
 Cd$-Kristallen. Z. 


3artschat, Optische und elektrische Eigenschaften v 
Physik 132, 224 (1953). 


i 


wurd 
Abb. 
desor 
a0 + 
l. 
klein 
2. 
als 6 
3. 
lel zu 
4. 
gestö 
gab. 
E 
zahl 
störu 
stark 
flexe 
D 
Krün 
stimr 
halb 
läßt 
Foku 
mung 
mung 


Für « 
der 7 

E 
übers 
zeich 
grenz 
Kriin 

A 
Gitte 
liche: 
sorte 
Abst: 
unten 
tarze 
mit 
dann 


10 
zu be 


| 
= 


Wolf. 
ird im 
S’, der 
let der 
ıkel a 
1 O der 
1 und 
einen 
ahlen- 
winkel 
und 
Da es 
t, wird 
Strah- 
einem 
ößeren 
diver- 
achher 
-hnung 
>benen 
‚ebene. 
kussie- 


re Wisse 


nfleck-| läßt sich bei kontinuierlicher Änderung der Orientierung immer ein neuer 
tgegenf Fokussierungspunkt finden. Unter Voraussetzung einer kreisförmigen Kriim- 
nicht}. mung kann man aus der Änderung Ad des Fokussierungsabstandes d den Krüm- 
sorien-| mungsradius R näherungsweise bestimmen: 
1g OF t2dad 
jeweil id 
ür ein Ad von 10 mm ergibt sich ein R von 200 m. Die Krümmungsradien 
der 7 gekrümmten CdS-Einkristalle lagen zwischen 10 und 150 m. 
Es ist anzunehmen, daß sich im Gitter 
10 mn} überschüssige Stufenversetzungen eines Vor- 
züchte' | zeichens ausbilden, die zu Kleinwinkelkorn- 
rdnung} grenzen aneinandergereiht die beobachtete 
en eim| Krümmung hervorrufen können. 
| zuein- Abb. 3 zeigt eine Projektion des CdS- 
ite def Gitters parallel zur c-Achse mit einer mög- 
ond die} lichen Stufenversetzung. Die beiden Atom- 
den jef sorten liegen in verschiedenen Ebenen im R 
breitenf Abstand c,/6 übereinander. Im rechten $ 
unteren Teilbild ist der Grundrißder Elemen- Abb. 3. Stufenversetzung im CdS- 
rgent. | tarzelle gestrichelt angedeutet. Siehtmandie itter; Projektion parallel zur c- 
CdS.) mit a bezeichnete Achse als a,-Achse an. 
1 der Elementarzelle; b Richtung : 
-istallen f dann kommt die Versetzung dadurch zu- des Burgersvektors 
ten Vi ) Zur Bestimmung der Desorientierungen ist jeweils die effektive Brennfleckbreite 


. 
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wurden mikroskopisch und verschiedentlich auch photometrisch bestimmt 1°). 
Abb. 2 zeigt Photometerkurven zweier verschieden breiter Reflexe mit Winkel- 


desorientierungen bzw. 
50 + 15 sec. 

Unter 30 Kristallen befanden sich: 

1. 13 Kristalle, deren Desorientierungen 
kleiner als 30 + 15 see sind, 

2, 22, deren Desorientierungen 
als 60 + 15 see sind, 

3. 7 Kristalle, die um eine Achse paral- 
lel zur c-Achse gekriimmt sind und 

4. ein Kristall, dessen Gitter sehr stark 
gestört war und völlig verzerrte Reflexe 
gab. 

Ein Zusammenhang zwischen der An- 
zahl der Rippen und dem Grad der Gitter- 
störung wurde nicht gefunden, auch sehr 
stark gerippte Kristalle können scharfe Re- 
flexe liefern. 


von 25+ 15sec 


kleiner 


Das Verfahren gestattet auch, die 
Krümmungsradien der Kristalle zu be- 
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Abb. 2. Photometerkurven der 
(0110)-Reflexe zweier verschieden 
stark desorientierter Kristalle 


stimmen. Während bei diskontinuierlicher Änderung der Orientierung inner- 
halb einer Netzebenenschar die Fokussierungseigenschaften verlorengehen. 


zu berücksichtigen. 
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stande, daß im oberen Teilbild eine (1010)-Ebene zusätzlich eingeführt ist. Der 
Burgersvektor hat in diesem Falle die Richtung der a,-Achse und sein Be- 
trag ist gerade 4,13 A. Aus dem bekannten Krümmungsradius R “3 dem 


Betrag des Burgersvektors ergibt sich eine Versetzungsdichte n = pdb. 
_ für die gefundenen Kriimmungsradien Versetzungsdichten von 1500 bis 2 2 5.000 
Versetzungen pro cm?. 

Nach Cullity und Julien!!) kann man mit der Guinier-Tennevin- 
schen Methode auch Torsionen im Kristall nachweisen, wenn man in den 
_ Primarstrahl, möglichst nahe der Strahlenquelle F, eine Sollerblende (SB) 
‚stellt. Dadurch wird der Strahlengang in einzelne Strahlungskeile zerlegt, 
die jeweils für sich einen Reflex erzeugen. Ist der Kristall ungestört, so liegen 
diese Reflexe alle an der gleichen Stelle und es tritt keine Änderung gegenüber 
de »m Aufbau ohne Sollerble nde (SB) auf. Liegt aber eine Torsion im Kristall vor, 
so wird einmal der Reflex aus seiner ursprünglichen Lage herausgedreht 
und zum anderen erzeugen die einzelnen Zellen der SB jede für sich ein Bild. 
Mit steigender Torsion wächst der Abstand zwischen den Einzelreflexen, und 
sowohl aus diesem Effekt als auch aus dem Winkel, um den der Reflex des 
tordierten Kristalles aus seiner ursprünglichen Lage herausgedreht wird, läßt 
sich der Torsionsgrad bestimmen. 

Abb. 4 (a, b. c) zeigt die künstlich hervorgerufene Torsion eines Kristalles 
von 0—4,8 min/mm Torsionsachse. Es handelt sich bei den Bildern um Reflexe 
der (0001)-Ebene. Die Torsion verteilt sich nicht gleichmäßig über den ganzen 
Kristall, wahrscheinlich infolge der nicht einheitlichen Dicke der zwischen 
benachbarten Rippen liegenden Abschnitte. 

In Abb. 5 (a, b) ist der (0001)-Reflex eines 0,2 mm dicken Kristalles ohne 
und mit SB aufgenommen. Das linke Ende entspricht der Aufklebestelle des 
Kristalles. Da nahezu bei allen Exemplaren eine solche Aufspaltung. wenn auch 


Abb. 4a Abb. 4b Abb. 4 


nach Cullity und Julien. Torsion von 0—4,8 min/mm Torsionsachse 


nicht immer so stark wie hier, auftrat, kann man daraus schlieBen, daB das 
Aufkleben das Kristallgitter oft betrachtlich verzerrt. Hier handelt es sich. 
wie Aufnahme b zeigt, um eine Torsion. Die Torsionsachse liegt in der (0001)- 
Ebene und ist parallel zur Bildebene. Die Torsion von oberem und unterem 
Kristallteil erfolgt nach entgegengesetzten Richtungen und betragt fiir den 
oberen Teil 2 min/mm und fiir den unteren Teil 0,9 min/mm Torsionsachse. 
Der Mittelteil des Reflexes deutet auf geringere Störungen hin, während die 


1), B.D. Cullityu.C. A. Julien, Effect of a Soller Slit on the Diffraction of X-Rays. 
J. appl. Phys. 24, 541 (1953). 
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rechte Seite wieder sehr starke Torsionen enthält, die beim Wachstum ent- 
standen sein müssen. Der Torsionsgrad läßt sich nur bestimmen. wenn man 
die Richtung der Torsionsachse kennt. Steht diese senkrecht auf der unter- 
suchten Netzebene, so tritt keine Reflexaufspaltung auf. Wie die der Abb. 5b 


entsprechende Aufnahme für die (0110)-Ebene zeigt, sind auch diese Ebenen 
tordiert. Das bedeutet, daß die Torsionsachse zwischen den beiden Ebenen 
liegt. Betrachtet man die wellenförmig verlaufenden Streifensysteme in 
Abb. 5b, so gewinnt man 
den Eindruck, als ob im 
rechten Kristallteil das 
Gitter einmal um volle 
360° tordiert ist. Mit 
Sicherheit ist das aber 
nicht festzustellen, da die 
genaue Lage der Torsions- 
achse nicht bekannt ist. 


Abb. 5. Im F-Punkt aufgenommener (0001)-Reflex eines 
Aufschlüsse über die Real- stark desorientierten Kristalles. 
struktur von Kristall- a) ohne b) mit Sollerblende nahe der Strahlenquelle 
blittchen erhält man, stehend 


wenn man den Reflex 
nicht im Fokussierungs- 
punkt, sondern méglichst 
dicht hinter dem Kristall 
auffängt. Bei einem 
Braggwinkel von 6 
haben sich nach 100 mm 
primarer und reflektierter 
Strahl schon weit genug Abb. 6a Abb. 6b 
getrennt, um eine solche Abb. 6. Der gleiche Reflex wie in Abb. 5 dicht hinter 
Aufnahme machen zu dem Kristall aufgefangen. 
nn Wises au a Dicht vor dem Kristall passierte die Strahlung die 
: Sollerblende mit zu den Netzebenen parallelen Spalten. 
einen Reflex, der die Um- b) Dicht vor dem Kristal! befand sich ein Sieb im Strah- 


risse des Kristalles hat, lengang 

und in dem sich die 

Fehler des Kristallinneren durch Schwärzungsänderungen bemerkbar machen. 
Wegen der geringen Konvergenz der Strahlung (Öffnungswinkel etwa 35 min). 
haben die Bilder Ähnlichkeit mit denen, die z. B. Lang!?) und Barth und 
Hosemann») mit der Parallelstrahlmethode erhalten. Da diese Methode im 
Gegensatz zu unserem Aufbau mit monochromatischer Strahlung arbeitet. 
ist sie jedoch wesentlich empfindlicher und gestattet bei nicht zu hohen 
Dichten den Nachweis einzelner Versetzungen. Unsere Bilder dienten ledig- 
lich dazu, bei schlecht fokussierenden Kristallen die im Fokussierungspunkt 
gemachten Aufnahmen zu deuten. 


12) A. R. Lang, Studies of Individual Dislocations . . . J. appl. Phys. 37, 1748 
(1959). 

13) H. Barth u. R. Hosemann, Anwendung der Parallelstrahlmethode im Durch- 
strahlungsfall zur Prüfung des Kristallinneren .... Z. Naturforschg. 13a, 792 (1958). 
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Die Abb. 6 (a. b) zeigt solche ,.Nahaufnahmen* des gleichen Kristalles 
wie Abb. 5. Dabei war in Abb. 6a die SB mit parallel zu den reflektierenden 
Netzebenen verlaufenden Spalten und in Abb. 6b ein feines Sieb kurz vor dem 
Kristall in den Primärstrahl gestellt worden. Im Reflex eines ungestörten 
Kristalles müßten sich die Muster ungestört abzeichnen. Hier deuten sie die 
Spannungen im Kristall an. Das rechte Reflexende der Abb. 6a scheint den 
schon aus Abb. 5b gezogenen Schluß zu bestätigen. daß der entsprechende 
Kristallbereich eine Torsion um volle 360° um eine nicht in der untersuchten 
Netzebene liegende Achse enthält. Die beiden im Reflex vertikal verlaufenden 
dunklen Streifen entsprechen je einer starken Rippe im Kristall. Der dreieck- 
förmige Punkt am linken Reflexende gibt an, wie weit sich die durch den 
eintrocknenden Lack nach entgegengesetzten Richtungen tordierten Bereiche 
ins Kristallinnere erstrecken. 

Ein möglicher Einfluß von elastischen Spannungen auf die Leitfähigkeit 
von CdS-Einkristallen sowohl bei optischer Grundgitter- und Ausläuferab- 
sorption als auch bei Röntgenanregung lag unterhalb der Meßgenauigkeit 
von 1%. 


Herrn Professor Dr. H. Simon danken wir für die Anregung zu dieser 
Arbeit und Herrn Dr. F. Eckart für anregende Diskussionen. 


Berlin. Physikalisch-Technisches Institut der Deutschen Akademie der 
Wissenschaften. Bereich elektronische Halbleiter. 


wh pot a4 3ei der Redaktion eingegangen am 26. Juli 1960. 
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Ober die Längenänderung 
von Ferromagnetikum-Träg erstoff-Gemischen 
im 


Von @. Vogler 


Mit 10 Abbildungen 


Inhaltsübersicht 


Es wird aus einer Extremalbetrachtung eine Beziehung für die Längen- 
änderungen, genannt Längsstreufeldstriktion, von Ferromagnetikum-Träger- 
stoff-Gemischen im Magnetfeld abgeleitet. Um die Richtigkeit dieser Abschät- 
zung nachzupriifen, wurden aus Carbonyleisen-Pulver und Gummilösung 
Proben hergestellt und deren Längenänderungen im Magnetfeld gemessen. 
Der Verlauf der an Zylindern gemessenen Kurven der Längenänderung in 
Abhängigkeit vom relativen Eisengehalt und für verschiedene Formfaktoren 
p stimmt qualitativ mit dem für Ellipsoide berechneten überein. Auf Grund 
der vorliegenden Untersuchungen ist es möglich, magnetische Werkstoffe 
herzustellen, die eine gewünschte Längenänderung unter Einwirkung eines 
bestimmten Magnetfeldes aufweisen. Diese Möglichkeit, die in diesem Umfang 
bei der Magnetostriktion homogener, kompakter Ferromagnetika nicht be- 
steht, dürfte zu mancher Anwendung Anlaß geben. 


Vor einigen Jahren veröffentli "ten Rust!), Rust und Pilz?) sowie 
Stauß und Sandoz?) Messungen der Volumen- und Längenänderungen von 
Ferromagnetikum-Trägerstoff-Gemischen im Magnetfeld. Während Rust 
und Mitarbeiter glaubten, die außerordentlich großen Volumenänderungen 
ihrer Ol-Ferromagnetikum-Aufschlimmungen auf Magnetostriktion homo- 
gener Ferromagnetika zurückführen zu können, erklärten Stauß und Sandoz 
die Längenänderung ihrer Probe mit der Drehung der in einer Gummimatrix 
eingebetteten Eisenflocken. 

Eigene Untersuchungen ®) ließen die abnormen Längen- und Volumenände- 
rungen der Gemische im Magnetfeld als Folge einer magnetischen Wechsel- 
wirkung zwischen den eingebetteten ferromagnetischen Mikroteilchen erkennen. 
die ihrerseits eine elastische Deformation des Trägerstoffes bewirkt. In einem 
Magnetfeld H, bildet sich um ein homogenes Ferromagnetikum ein Streufeld 


1) H. Rust, Z. angew. Physik 3, 9 (1951). 

2) H. Rust u. P. Pilz, Z. angew. Physik 3, 379 (1951). 

3) H.E. Stauß u. G. Sandoz, J. Appl. Phys. 25, 267 (1954). 
4) (1. Vogler, Z. physik. Chem. 205, 360 (1956). 
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aus, das von der Gestalt des Ferromagnetikums abhängig und dem äußeren 
Feld H überlagert ist. Bei den Mischkörpern ist diesen beiden Feldern auBer- 
dem noch das stark inhomogene Streufeld der Mikroteilchen überlagert. 
Beide Streufelder — Makrofeld und Mikrofeld — sind mit einer Streufeldenergie 
verknüpft). 

Die Mischkörper deformieren sich in einem äußeren Feld so, daß die Summe 
aus magnetischer und elastischer Energie ein Minimum wird. Dabei sind in 
Abhängigkeit vom Dimensionsverhältnis sowohl negative als auch positive 
Längenänderungen möglich. In erster Näherung ist anzunehmen, daß der 
Längenänderung eine gleichsinnige Volumenänderung entspricht, wenn die 
Poissonsche Zahl kleiner als 0,5 ist. 

Zur Abschätzung der Längenänderungen der Ferromagnetikum-Träger- 
stoff-Gemische im Magnetfeld soll eine Modellsubstanz benutzt werden, bei 
der kugelförmige, ferromagnetische Teilchen in einer elastischen Trägersub- 
stanz statistisch verteilt eingebettet sind. Damit ist der unübersichtliche 
Effekt der Drehung der Teilchen im Magnetfeld ausgeschlossen, der bei läng- 
lichen Teilchen berücksichtigt werden müßte. 


a Theorie der Längsstreufeldstriktion 


Der Berechi ial der Längenänderung soll ein ellipsoidförmiger Körper 
zugrunde gelegt werden. In einem homogenen Feld deformiert sich ein Körper 
so, daß die Gesamtenergie (Feldenergie, äußere Streufeldenergie. innere Streu- 
feldenergie, elastische Energie) zum Minimum wird. 


E,=E£,+£,4+ + = Min 
Bg = — Hedy Va +4 No Jj Vo Mi + 
E, — Gesamtenergie H äußeres Feld 
, — Feldenergie J, — pauschale Magnetisierung 
E, — äußere Streufeldenergie N „— äußerer Entmagnetisierungs- 
faktor 
E, — innere Streufeldenergie N. — innerer Entmagnetisierungs- 
faktor 
Ezı — elastische Energie N „ — Entmagnetisierungsfaktor der 
Kugel 
la — Längenänderung 6 — Elastizitätsmodul 
A — Lange der Probe ohne Feld uu — effektive Querzahl 
by — Durchmesser der Probe ohne 7 
Feld 
p=; Formfaktor ohne Feld 
’o 
— Gesamtvolumen 
v — relatives Ferromagnetikum- 


volumen 
Indem man den Ausdruck fiir die Gesamtenergie nach a ableitet und Null 
setzt, erhalt man das Minimum der Gesamtenergie und durch Auflésen nach 
la/a, die gesuchte relative Längenänderung des Körpers im Magnetfeld. Bei 


der Ableitung des Ausdruckes der Gesamtenergie wird angenommen, daß sich 


5) G. Vogler, Naturwiss. 13, 423 (1959). E ‘+ 
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in erster Näherung die statistische Verteilung bei einer Längenänderung des 
Körpers nicht ändert. 


da föp EN, I. oJ, N,J, 
A, ~ éa op [- 2 Ji Tan, (- — Nady +H, 


‘ty v 
+ cv v Na Jot H,) at ov ] 
eV, N,J2 H,J,)\ 


Die Schreibweise der partiellen Ableitungen ( - ‘) und ( =), die korrekt 
cv cv 


Ov 
lichkeit beibehalten werden. Außerdem wurden höhere Potenzen von Aa ver- 
nachlässigt. 

Mit 


const const 


ad» 
und ( —3 ), lauten müßten, soll aus Gründen der Übersicht- 
cv Nal 
La 


ao 


_9 
ca 1-24), da (I 2 


ao 


erhält man für die relative Längenänderung von Ferromagnetikum-Trager- 
stoff-Gemischen : 


la oN, 2 (_ J, ON; N; 
Ip ( éN, N.— J, ] (1) 


2vaJ,(N; H\ 2N; yy , 2H 

2a) | ev ( N, 2.) v J ( ev 

In der Sättigung (J, = J,,) ist N x 


[V; = N, (1-—v)] und j (- <1. 


Damit vereinfacht sich die Formel für die relative Längenänderung und 
lautet für die Sättigung: 


la J2, 
= 


2. Probenherstellung und Messung 

Die Messung der Längenänderungen im Magnetfeld erfolgte an zylindri- 
schen Mischkörpern aus Gummi und Eisen, die dem bei der Abschätzung der 
magnetischen Streufeldstriktion benutzten Kugelmodell weitestgehend ent- 
sprachen. Abgemessene Volumina Gummilösung und Carbonyleisen-Pulver C 
(Leuna-Werke) wurden zusammengeschüttet, gut verrührt und, nachdem ein 
großer Teil des Lösungsmittels verdunstet war, mit Hilfe eines Walzwerkes 
geknetet und zu homogenen, blasenfreien, zylindrischen Proben verarbeitet. 
Brauchbare Körper können nach den bisherigen Erfahrungen nur mit Eisen- 
gehalten bis zu 50% hergestellt werden. Der relative Eisengehalt wurde durch 
Wiegen des Carbonyleisen-Pulvers .vor dem Mischen und durch Bestimmung 
des endgültigen Probenvolumens ermittelt. 

Von den so hergestellten Proben mit den relativen Eisengehalten von 9%. 
20%, 31%, 40%, 50%. die einen Durchmesser von etwa 0,6 cm hatten, wurden 
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Proben mit 6 cm Länge abgeschnitten, gemessen, nach der Messung verkürzt äuß 
und wieder gemessen. Auf diese Art konnten Meßkörper mit den Formfak- der 
toren p = 10, p = 5, p = 3 und p = 2 hergestellt und ihre Striktionskurven wur 
ermittelt werden. 

Zur Messung der Streufeldstriktion stand eine im Institut fiir Magnetische 
Werkstoffe entwickelte Magnetostriktionsapparatur zur Verfiigung®), die auf 
einem Meßverfahren beruht, bei dem die Längenänderung der Probe auf die 
bewegliche Platte eines Kondensators übertragen wird, der sich in einer Brük- 
kenschaltung befindet. 

Die Gummiproben lassen sich wegen ihrer Konsistenz und der geringen 
Kräfte, die für die Bewegung des Meßsystems zur Verfügung stehen, nur 
schlecht messen, so daß die Meßfehler, verglichen mit denen bei Messungen 
an Metallen, sehr groß sind. 


3. Bestimmung der für die theoretische Abschätzung benötigten Größen 


Zur Prüfung der Theorie wird die Beziehung (1) benutzt. Diese enthält bis 
auf die Anpassungskonstante w, die die Bedeutung 2iner effektiven Querzahl 
hat, Größen, die aus dem Experiment oder aus Tabellen entnommen werden 
müssen. Die Messung der Größen J, und N, in Abhängigkeit von H, und v, 
über die im folgenden berichtet werden soll, ist deshalb nötig, weil der funk- 
tionelle Zusammenhang dieser Größen nicht bekannt ist. Obwohl die Mes- t 


4 
N, 


7 200 400 600 800 1000 1200 1400 200 400 600 800 1000 1200 1400 
Abb. 1. Kurven der pauschalen Ma- Abb. 2. Scherungskurven der pau- 
gnetisierung für verschiedene relative schalen Magnetisierung einer Probe 
Eisengehalte v(V, = 0) mit 50°, Eisen in Abhängigkeit vom 
Feld. Parameter ist der äußere Ent- 
magnetisierungsfaktor N, 


sungen der Längenänderung Ja/a, an zylinderförmigen Proben vorgenommen 
wurde, soll die Berechnung der Streufeldstriktion für ellipsoidförmige Proben 
durchgeführt werden. Es ist begrifflich sauberer, wenn man, ebenso wie bei 
der Theorie, die Berechnung für homogen magnetisierte Körper durchführt s 
und nicht zum Teil inhomogen magnetisierte Zylinder zugrunde legt, deren 
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6) G. Vogler, Nachrichtentechnik 5, 310 (1955). 
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äußerer Entmagnetisierungsfaktor in verwickelter Weise von der Form und 
der Suszeptibilität der Probe abhängt. Die pauschale Magnetisierung (Abb. 1) 
wurde mittels einer ballistischen Meßanordnung ermittelt. Aus den für ver- 


600 


400 


200 2007 
4 iL i 
# or 02 03 % 
Abb. 3. Pauschale Magnetisierung Abb. 4. Pauschale Magnetisierung — 
einer Probe mit 50%, Eisen als Funk- als Funktion des relativen Eisenge- 
tion des äußeren Entmagnetisierungs- haltes. Parameter ist der äußere 
faktors. Parameter ist die Feld- Entmagnetisierungsfaktor 


stärke H, 


schiedene N, gescherten Kommutierungskurven J,(H,) (Abb. 2) wurden für 


asd 
jedes v die Kurvenscharen J,(N,) (Abb. 3) und daraus durch graphische 
oJ » ON, 
Differentation (; 7) entnommen. „” wurde aus den bekannten For- 
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05 > 
07 02 03 04 05% 
Abb. 5. Pauschaler innerer Entmagneti- Abb. 6. Pauschaler innerer Ent- 
sierungsfaktor als Funktion des relativen magnetisierungsfaktor bei H, 7 
Eisengehaltes. Parameter ist die Feld- 1400 Oe als Funktion des relativen 
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> 


mein für Ellipsoide errechnet und (% ) aus den J, (v)-Kurven (Abb. 4) 


const 

durch graphische Differentation ermittelt. Uber die Messung (Abb. 5) und 
Berechnung des inneren Entmagnetisierungsfaktors ist an anderer Stelle 
bereits berichtet worden 5)’). Die Werte für ‘) und (<4), wur- 

den aus den Diagrammen Abb. 6 und Abb. Tentnommen. Der Elastizitatsmodul 
ist in der Magnetostriktionsapparatur mittels einer Zusatzeinrichtung gemessen 


N und der Formfaktor p aus den geometrischen 
N: u Abmessungen berechnet und aus Tabellen das 

/ 0.09 N „ entnommen worden. 
2 a Nachdem durch Experimente sicher- 
gestellt worden war, daB sich bei einer Ver- 


™ 


m 040 kürzung der Probe die Querdimensionen 
0,50 vergrößern, also ähnliche Verhältnisse wie 
in der Mechanik vorliegen, wurde eine effek- 
tive Querzahl u = == eingeführt. Bei der 
magnetischen Streufeldstriktion wird die Län- 
genänderung der Probe nicht durch rein 


mechanische, sondern magnetische Kräfte 
Ip [6] hervorgerufen. Da die verformenden Kräfte 
nicht von außen angreifen, sondern im Inne- 
200 «00 600 800 ren der Probe wirken und recht unterschied- 
Abb. 7. Innerer Entmagnetisie-  ]Jiche Richtungen haben, soll hier von einer 
rungsfaktor als Funktion der pau- ; 
schalen Magnetisierung. Parameter effektiven Querzahl gesprochen werden. u 
ist der relative Eisengehalt wurde als einzige Größe benutzt, um die 
Abschätzung dem Experiment anzupassen. 
Dabei wurde es für jedes Mischungsverhältnis nur einmal passend gewählt und 
dann unabhängig von der Form des Körpers konstant gehalten. Aus der 
Anpassung der Abschätzung an die Messungen ergaben sich folgende u-Werte: 
vr = 0.09 u = 0,5 v = 0,40 u = 0,25 
v = 0,20 u 0,45 v = 0,50 u = 0,20 
v = 0,31 u = 0,3 
Die Abnahme des ~ mit wachsendem v ist verständlich, da bekanntlich die 
Poissonsche Zahl bei Gummi ohne Füllstoff, je nach dem Vulkanisations- 
grad, zwischen 0,5 für weichen und 0,2 für harten Gummi liegt. Die Zugabe 
von Eisenpulver hat auf den Gummi eine ähnliche Wirkung wie die stärkere 
Vulkanisation. Je fester und härter der Gummi wird, um so mehr nimmt die 
Poissonsche Zahl ab. Die einzige Anpassungskonstante der Abschätzung 
weist also einen Verlauf auf, der durch Messungen an eisenfreien Gummi 
bestätigt werden kann?). 


| 
I 


4. Vergleich und Diskussion der gemessenen und berechneten Längenänderungen 
In Abb. 8 ist für Zylinder mit kugelförmigen Mikroteilchen die relative 
Längenänderung Aa/a, für ein Feld von 1400 Oe über dem relativen Eisen- 


7) G. Vogler, Ann. Physik (7) 7, 268 (1961). 
8) J.D’ Ans u. E. Lax, Taschenbuch f. Chem. u. Phys., Springer-Verlag 
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volumen aufgetragen. Parameter ist der Formfaktor p der Mischkörper. Die 
Punkte für den gleichen Formfaktor p sind durch Geraden untereinander ver- 
bunden. Die für p = 10 (überwiegend innere Streufeldenergie) vorhandene 
Kontraktion wird mit abnehmendem p (zunehmende äußere Streufeldenergie) 
immer kleiner und geht schließlich in eine Dilatation über. In Abb. 9 sind 
für ellipsoidförmi- 
ge Proben mit | 
statistisch verteil- ——70% 
ten Mikroteilchen 
die mit der Be- 

ziehung (1) be- | 
rechneten relati- 800 7 
ven Längenände- 
rungen aufgetra- 6007 
gen. Beim Ver- 
gleich der gemes- 
senen und berech- 
neten Kurven 
stellt man fest, 
daß die Änderung 
der Streufeldstrik- 
tion in Abhängig- 
keit vom Form- 
faktor p und vom 


relativen Eisen- -400} 
gehalt v gut wie- . 
dergegeben wird, 01 02 03 04 05 01 02 03 04 05 


daß jedoch die Abb. 8. An Zylindern gemes- | Abb. 9. Berechnete relative 
sene relative Längenänderung Längenänderung für Ellipsoide 
p bei 1400 Oe in Abhängigkeit bei 1400 Oe in Abhängigkeit 
ches p nur 1M vomrelativen Eisenvolumenv. vom relativen Eisenvolumen v. 
Falle p= 2 und Parameterist der Formfaktorp Parameter ist der Formfaktor p 


p=53 innerhalb 

der Fehlerbereiche übereinstimmen. Besonders im Gebiet negativer Längen- 
änderungen sind die Abweichungen der berechneten von den gemessenen 
Werten groß. Sie finden ihre Erklärung in der Hauptsache darin, daß die 
Rechnung für Ellipsoide und die Messung an Zylindern durchgeführt wurde. 
Bekanntlich hat das Ellipsoid bei gleichem p einen größeren Entmagneti- 
sierungsfaktor als der Zylinder. Dieser Unterschied nimmt mit wachsendem 
pzu. Aus dieser Tatsache erklärt sich auch die für p= 2 und p = 3 gute 
und für p = 5 und p = 10 schlechte Übereinstimmung zwischen den gemes- 
senen und berechneten Werten. Deshalb vollzieht sich der Nulldurchgang 
der Streufeldstriktion beim Ellipsoid auch eher als beim Zylinder. So ist es 
zu erklären, daß bei den Zylindern mit p = 5 gerade der Nulldurchgang statt- 
findet, während die Ellipsoide mit p = 5 nur positive Streufeldstriktion auf- 
weisen. Die Abweichungen zwischen den berechneten und den gemessenen 
Werten der Längsstreufeldstriktion der Probe mit dem relativen Mikroteilchen- 
volumen von v = 0,09 können zu einem Teil auf die ‚‚negative Oberflächen- 


Kurven fiir glei- 


effektive Querzahl von u = 0,5 haben, beispielsweise bei der Probe mit 9%, 
26b 


1 i 
‘ 
7 
“oe 
fe’ 
a = 
spann ) dieser Prob ickgeführt werden. Bei Meßkörpern, die ein u 
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Eisenpulver, kann die negative Oberflächenspannung zu einer Vergrößerung 
_ der positiven Streufeldstriktion beitragen. 
In Abb. 10 ist der berechnete Verlauf der Längsstreufeldstriktion in der 
Sättigung für dieselben Werte von p und v wiein Abb. 9 aufgetragen. Aus Ver- 
suchen im Weißschen Magneten ergab sich, daß die Sättigung der Proben 
erst bei Feldern zwischen 10000 und 13000 Oe 
4 erreicht wird. Da fiir die Messung der Streufeld- 
’ striktion sehr homogene Felder benötigt werden, 
kam eine Messung der Längenänderung im Weiß. 
schen Magneten nicht in Frage. Eine Spule, die 
es gestattet, bei so hohen Feldern zu arbeiten, 
stand zur Zeit der Untersuchungen nicht zur 
Verfügung. Die Richtigkeit der Abschätzung 
muß daher zu einem späteren Termin nachgeprüft 
werden. 

Beim Vergleich der Abb. 9 und 10, in die 
die theoretischen Werte der Längsstreufeldstrik- 
tion für 1400 Oe und die Sättigung eingetragen 
sind, fällt auf, daß sich das Maximum der posi- 


4a 
4000 


10 


5000 


2000 


tiven Streufeldstriktion von v = 0,31 nach 
v = 0,2 verschoben hat und daß sich beispiels- 
p=0 weise für v = 0,31, p = 5 und v = 0,4, p=5 bei 


1400 Oe positive und in der Sättigung negative 
Längsstreufeldstriktionswerte ergeben. Über dem 

v äußeren Feld aufgetragen, erhält man bei Körpern 

dieser Form und Zusammensetzung also Kurven 

01 02 03 04. 05 der relativen Längenänderung, die den Magneto- 


Abb. 10. Berechnete relative striktionskurven des Eisens sehr ähnlich sind 

soide in der Sättigung in Ab- 
hängigkeit En Diese Abhängigkeit der Streufeldstriktion von 


Eisenvolumen. Parameter ist der Feldstärke H, kommt dadurch zustande, daß 
der Formfaktor p sich mit wachsender Feldstärke die N,-Kurven 
(Abb. 5) immer mehr der Sättigungsgeraden 
nähern. Mit dieser Näherung ist bei allen Proben eine Vergrößerung des inne- 

ren Entmagnetisierungsfaktors und, außer bei den Proben mit einem v < 0,2, 


_ eine merkliche Zunahme der Größe (5) verbunden. Während sich in For- 


mel (2) die positiven Glieder kaum ändern, nimmt bei Proben mit einer kleinen 
effektiven Querzahl u das negative Glied, das (1 — 2) als Faktor hat, stark 
zu und bewirkt dadurch eine Behinderung der Dilatation und zum Teil sogar 
eine Umkehr des Vorzeichens der Längsstreufeldstriktion. 


5. Einfluß des Materials auf die Längenänderung 

In den vorhergehenden Abschnitten war untersucht worden, wie sich die 
 Langsstreufeldstriktion mit dem Formfaktor p und dem relativen Eisengehalt 
ändert. Dabei wurde festgestellt, daß in dem untersuchten Bereich die Über- 
_ einstimmung zwischen der rechnerischen Abschätzung und den Messungen 
_ zufriedenstellend ist, wenn man berücksichtigt, daß nur eine grobe Näherung 
bei der Theorie Verwendung fand. Dabei blieb die Frage nach dem Einfluß 
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3erung | des verwendeten ferromagnetischen Materials unbeantwortet, obwohl zwei 
geometrische Größen — die Form der ferromagnetischen Teilchen und ihre 
in der | Richtungsverteilung — und eine magnetische Größe — die Sattigungsmagneti- 
is Ver. | sierung — einen wesentlichen Einfluß auf die Deformation des Meßkörpers 
-roben | haben. 
)00 Oe Vergleicht man Magnetostriktionskurven mit Längsstreufeldstriktions- 
ufeld. | kurven, so fällt auf, daß zwar die Längsstreufeldstriktion in der Sättigung 
erden, | wesentlich größer ist als die bekannten Sättigungsmagnetostriktionswerte, daß 
Neiß- | aber sehr große Felder zur Sättigung benötigt werden. Das ist in erster Linie 
le, die | darauf zurückzuführen, daß die Kurven der pauschalen Magnetisierung (Abb. 1) 
beiten, | bei kugelförmigen Teilchen wegen des großen inneren Entmagnetisierungs- 
ıt zur | faktors zu stark geschert sind. Wenn es gelänge, elliptische Teilchen ausgerich- 
itzung | tet anzuordnen, dann würde man den inneren Entmagnetisierungsfaktor und 
eprüft | damit die innere Streufeldenergie herabsetzen und bei kleineren Feldern die 
Sättigung erreichen. 
in die Da in die Beziehungen für die Längsstreufeldstriktion die pauschale 
istrik- | Magnetisierung quadratisch eingeht, könnte man bei Verwendung von Mikro- 
tragen | teilchen aus einer geeigneten Eisen-Kobalt-Legierung an Stelle von Carbonyl- 
" posi- | eisen-Pulver die Sattigungsstreufeldstriktion um 25% steigern. 

nach Da der für die Versuche verwendete Gummilösungsgummi für eine An- 
spiels- | wendung der magnetischen Streufeldstriktion nicht geeignet ist, wurde eine 
- 5 bei | Anzahl Kunststoffe untersucht, die genügend elastisch sind und denen man 
gative | einen größeren Prozentsatz Eisen-Pulver als Füllstoff zugeben kann. Die 
r dem | Untersuchungen ergaben, daß man Buna S 4 bis zu 50%, Eisen zugeben kann 
örpern | und daß sich durch verschieden starke Vulkanisation der E-Modul in weiten 
‘urven | Grenzen ändern läßt. Auf diese Weise ist es möglich, in Abhängigkeit von der 
gneto- | Probengestalt, dem relativen Eisengehalt, dem Vulkanisationsgrad und dem 
sind | äußeren Feld jede beliebige relative Längenänderung zwischen — 400: 10% und 
+ 1000. 10-8 hervorzurufen. Abschließend sei noch betont, daß es durchaus 
»n von | möglich ist, daß sich dieser Werkstoff besonders gut als Unterwasserschall- 
je, daß | geber oder Empfänger eignet. 


‚urven 
eraden 6. Streufeldstriktion bei Ferriten und Dielektrika 
; inne- Nach den bisherigen Betrachtungen erhebt sich die Frage, ob die Streufeld- 


< 0,2, | striktion einen Einfluß auf die gemessenen Werte der Magnetostriktion hat. 
n For- | Die äußere Streufeldenergie macht sich bekanntlich?) bei ungünstig geformten 
. Proben so stark bemerkbar, daß sie zu Meßfehlern Anlaß geben kann. Diese 
leinen | Fehler lassen sich innerhalb der Meßgenauigkeit vermeiden, wenn man zur 
stark Messung der Magnetostriktion nur lange, dünne Proben verwendet. Die innere 
‚Sogar | Streufeldstriktion kann bei metallischen Proben nur dann einen merklichen 
Einfluß auf die Magnetostriktionsmessungen haben, wenn es sich um hetero- 
genes Material handelt, das Phasen unterschiedlicher Magnetisierung enthält 
und dessen Elastizitätsmodul klein ist. 
ch die Bei den Ferriten, die unter Umständen zahlreiche Hohlräume enthalten, 
gehalt besteht die Möglichkeit, daß die Streufeldstriktion die Magnetostriktions- 
Über- messungen verfälscht. Da die äußere Streufeldenergie durch günstige Form- 
sungen | gebung stark verringert werden kann, wurde nur der Einfluß der inneren 
1erung 
'influß 9) Becker-Döring, Ferromagnetismus, Springer-Verlag. 
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> Streufeldenergie abgeschätzt. Dabei ergab sich, daß unter extrem ungünstigen 


Verhältnissen — viele Löcher, große Magnetisierung, kleiner Elastizitätsmo- 
dul — maximale relative Längenänderungen von — 1- 107 zu erwarten sind. 
Es ist anzunehmen, daß im allgemeinen die durch die Streufeldstriktion er- 
zeugte Längenänderung mindestens eine Zehnerpotenz unter der berechneten 
oberen Schranke liegt, so daß dieser Effekt bei der Magnetostriktionsmessung 

an Ferriten nicht berücksichtigt zu werden braucht. 
Es ist selbstverständlich, daß es der magnetischen Streufeldstriktion ana- 
_ loge Erscheinungen bei den Dielektrika im elektrischen Feld gibt (dielektrische 
_ Streufeldstriktion). Allerdings kommt hier als wesentlicher Unterschied hinzu, 
daß die Dielektrizitätskonstante des Trägerstoffes, die immer kleiner als die 
der Teilchen sein soll, berücksichtigt werden muß und man bei der Elektri- 
sierung — Ferroelektrika ausgenommen — den Zustand der Sättigung nicht 
kennt. Mit wachsender Dielektrizitätskonstante des Trägerstoffes nimmt die 
innere Streufeldenerzie ab, um — isotropes Dielektrikum vorausgesetzt — den 
Wert Null anzunehmen, wenn beide Dielektrizitätskonstanten übereinstim- 
_ men. Je größer der Unterschied der Dielektrizitätskonstanten der beiden 


stoffes gegenüber der der eingebetteten Teilchen vernachlässigen. Besonders 
gute Übereinstimmung zwischen dielektrischer und magnetischer Streufeld- 
striktion müßte sich mit Proben erzielen lassen. bei denen ferroelektrische 
Teilchen in einem Kunststoff mit sehr kleiner Dielektrizitätskonstante einge- 
 bettet sind. Abschließend soll noch erwähnt werden, daß es möglich sein 
müßte, Proben herzustellen, die sowohl magnetische als auch dielektrische 
_ Streufeldstriktion aufweisen. Bekanntlich findet man bei den meisten Ferriten 
eine hohe Dielektrizitätskonstante mit guten magnetischen Eigenschaften 
vereint, so daß Ferrit-Pulver-Trägerstoff-Gemische im magnetischen und 

elektsischen Feld ihre Form ändern. 


‘ Herrn Prof. Dr. W. Schiitz habe ich fiir die Unterstiitzung bei der Durch- 
führung der Untersuchungen ganz besonders zu danken. 


Jena, Institut für Magnetische Werkstoffe der Deutschen Akademie der 
zu Berlin. 
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Berechnung der 


entarteter anharmonischer Oszillatoren 
mit Hilfe der Mellinschen Integraltransformation = 


Inhaltsiibersicht 


Es wird ein neues Verfahren zur Berechnung der Zustandssumme eines 
mechanisch entarteten anharmonischen Oszillators angegeben, das auf der Ver- 
wendung der Mellinschen Integraltransformation beruht und auch fiir die 
Summation anderer unendlicher Reihen, insbesondere solcher mit Fakultäten, 
benutzbar ist. Als Anwendungsbeispiel wird die Zustandssumme eines Systems 
entarteter anharmonischer Oszillatoren in &-Schreibweise (w;: Frequenz 

[em-!] für unendlich kleine Amplitude) berechnet. rae 


da 


Bei der Berechnung der Zustandssumme eines mechanisch entarteten 
anharmonischen Oszillators mit dem Entartungsparameter d treten bekannt- 
lich Summen folgender Form auf: 


1777 


Hierin bedeutet n die Schwingungsquantenzahl, und es ist k > 0 und ganz- 


zahlig; gp = (" = stellt das statische Gewicht des Zu- 
hco 


standes mit der Schwingungsquantenzahl n dar; weiter ist u = LT mit 


o = Frequenz (in cm~) des Oszillators für unendlich kleine Amplitude. Für 
derartige Summen existiert ein allgemeines Berechnungsverfahren von Wool- 
ley!), bei dem durch Differentiationen und verwickeltere algebraische Um- 


formungen und Substitutionen von der Summe x" ausgehend die gewünsch- 
> 
n 


ten Summenformeln erhalten werden. Es wird hier ein anderes Verfahren 

gezeigt, das die Mellinsche Integraltransformation, eine Abart der bekann- 

teren Laplace-Transformation, verwendet. Dieses Verfahren kann auch auf 

andere Reihen als (1) angewandt werden und eignet sich besonders zur Sum- 

mation von Reihen, deren Glieder Fakultäten enthalten. Ein entsprechendes 

Verfahren zur Summation anderer Reihen unter Verwendung der Laplace- 
Transformation ist von Wheelon angegeben worden?). 

') H. W. Woolley, Res. Nat. Bur. Stand. 56, 105 (1956). 

2) A.D. Wheelon, J Phys. 25, 113 (1954). 
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Die Mellin-Transformierte M{f(x), n} = f(n) einer Funktion f(x) ist defi. 
niert als?) 

fin) = M{f (x); n} = f f(x) a®-1dz, (2) 
0 


vgl. Sneddon‘). Andererseits gilt fiir die Gammafunktion die Definitions. 


4 0 = 
Es liegt daher nahe, zur Berechnung von Summen wie er 
— vgl. (1) — die Mellin-Transformation als Hilfsmittel heranzuziehen. (Die 


Laplace-Transformation selbst scheint fiir derartige Aufgaben nicht geeignet 
zu sein.) 


Der Grundgedanke des Summationsverfahrens ist der folgende. Man 
schreibt die zu berechnende Summe in der Form 


- 2, a, f(n), () 


wobei /(n) die Mellin-Transformierte einer Funktion f(x) und a, ein noch 
von n gesetzmäßig abhängender Koeffizient ist. Unter der Voraussetzung der 
Vertauschbarkeit von Summation und Integration wird 


n=0 n=0 fi) 

Der allgemeine Summand der zu berechnenden Summe S ist möglichst ge- 
schickt in die beiden Faktoren a, und f(r) aufzuspalten, so daß sowohl die 
Summe a, x"! wie das Integral f f(x)... dx geschlossen ausgewertet werden 
können. 


Zur Berechnung der Summe (4) gibt es verschiedene Möglichkeiten. Ein 
einfaches Verfahren beruht auf der Verwendung der beiden Beziehungen‘) 


von denen die erste nach (3) die Definition der Funktion J’(n + d) darstellt 
und die zweite leicht durch wiederholte partielle Integration bewiesen werden 


3) Eine Verwechslung des Differentialsymbols d mit dem Entartungsparameter 4 
ist wohl ausgeschlossen. 
4) I. N. Sneddon, Functional Analysis, in Handbuch der Physik, Bd. II, 1955. 
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kann. Damit ist 
n=0 


mit 
-un 


éx J 
Nach (6) ergibt sich 


Die Summe wird mit x” = e"!n® zu 


m+d-Di 5 
n! — 


(10) 


und 


© (‚Inz-u)" 
(13) 
i=. 


und es entsteht die endgiiltige Beziehung 


n ‘= 0 n! 


Die hier auftretenden Integrale lassen sich sofort mittels 


co 
f ge = (15) 


geschlossen auswerten. Man erhält z.B. firk=1 


co 


Ein anderes Verfahren zur Berechnung der Summe (4) bzw. (9) läuft über den 
Ansatz 


k 
a) on! 


f(n) = (n +d—1)! = M {ate-*; n}. 


n* 
A, =e lim (5 


(A-u)n 
f(r) = S lim (5) (n+d-—1)! 
el4- u+lnz)n (19) 
= lim (5 —— dx 
a 0 n= 0 wi 
= -tyd-1 (e x) 


| 
ay (n) L (9) | 
= 
. Mar 
6 = 
r 
6 
=, 
- 
= 
_ 
Man erhält 
4 
i= 
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und es ergibt sich die endgültige Formel ris 


In.(19) muß die Vertauschbarkeit von Summation und Grenzwertbildung und 
danach die Vertauschbarkeit von Summation und Integration vorausgesetzt 
werden. 


IE Das eben gezeigte Verfahren mit der Einführung des Parameters A und 
FE dem Grenzübergang 4 — 0 ist bei vielen Summationsaufgaben von Nutzen. 
En Bei der Summation der Reihe (4) ist die Verwendung eines Parameters A 
j wegen des Vorhandenseins der Exponentialfunktion überflüssig, wenn die 


Summe für k = 0, 


= (n +d—1)! eun— - (d — 1)! (21) 
als bekannt vorausgesetzt wird. Es ist nämlich 
= (— (=) (22) 
. (n+ 4-1)! nkeun — (— 1)* 2) (n+d—1)! 
(23) 


_ Die endgültige Beziehung lautet jetzt 


Für beliebiges d und k = 0,...,3 ergeben sich folgende Formeln: bg 4 
(25) 
x (" +d- ') Ben 
* 


WwW 
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- _ Als Anwendungsbeispiel werde zum Schluß die Zustandssumme eines Systems = 
entarteter anharmonischer Oszillatoren berechnet. Durch ein solches System 
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wird z. B. das Schwingungsverhalten der Atome eines mehratomigen Molekiils 
beschrieben. Bei anharmonischen Schwingungen wird in der Molekülspektro- 
skopie aus praktischen Gründen meist die v-Schreibweise verwendet; das Sym- 
bol »; bedeutet eine beobachtete Frequenz (in cem-!), also die Differenz zweier 
Termwerte (vgl. Herzberg®)). Für theoretische Betrachtungen ist dagegen 
die w-Schreibweise nützlich; die Größe w, ist die entsprechende Frequenz 
(in em-!) für unendlich kleine Amplitude, also für harmonische Schwingungen. 
Die Schwingungsenergie des betrachteten Systems ohne den Anteil der Null- 
punktsenergie sei mit x,.: Anharmonizitätskonstante (in cem-!)®) 


2; +22 N; +> (4; m + d, n) |}; (26) 


der Anteil des Drehimpulses der bei mechanisch entarteten Schwingungen 
möglichen Rotationsbewegungen der Atome ist der Einfachheit halber wegge- 
lassen. 

Die Zustandssumme des Systems ohne den Anteil der Nullpunktsenergie 
wird 


P 
N, Na i N; 
Mit uv=hcw,kT und y.=hcx,/k T wird 


E 1 
= un; + = + d;,) + >) Yır N. (d; + 4, |. 


Einsetzen von (28) in (27) und Verwendung der Abkürzung g,; = re +4~ ‘) 
ergibt 


k(<i) 


Q II Io e 


I 


= [I + dim, 
Yix + 4 + d,n,)]\ {| 
>» Yi + 4 (d,n, +d,n, di}. 


Die Summation ist nun mit den Beziehungen (25) durchführbar ns 4 ergibt 
nach einiger Rechnung 


5) G. Herzberg, Infrared and Raman Spectra of Polyatomic Molecules, 1945. 
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Das Produkt stellt die Zustandssumme der harmonischen Schwingungen dar, 


_ der Faktor in geschweiften Klammern ist die Anharmonizitätskorrektur. 


Der Verfasser möchte auch hier Herrn Prof. Mühlenpfordt für sein 
stetes Interesse und großzügiges Entgegenkommen danken. 


Leipzig, Arbeissstelle für Statistische Physik der Deutschen Akademie 
der Wissenschaften zu Berlin. 
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Ein neuer Ansatz 
für die Behandlung gasdynamischer Probleme 
bei starken Abweichungen 
vom Thermodynamischen Gleichgewicht 


Von Fritz Weitzsch 


Inhaltsübersicht 


Es wird ein neuer Ansatz zur Gewinnung von Näherungslösungen der 
Boltzmann-Gleichung vorgeschlagen, der gegenüber den bisherigen Ansätzen 
eine Reihe von Vorteilen bietet. Die im Ansatz verwendete Geschwindigkeits- 
verteilung ermöglicht in der Hauptsache eine bessere Anpassung an das je- 
weils vorliegende gasdynamische Problem, bzw. an die jeweils zu erwartenden 
Abweichungen der Geschwindigkeitsverteilung von der Gleichgewichtsver- 
teilung. 

Die Brauchbarkeit des Ansatzes wird am Beispiel des stationären Ver- 
diehtungsstoßes geprüft. Bei annehmbarem Aufwand werden modelltheore- 
tisch beliebig starke Stöße untersucht. Die Lösungen können als nächst- 
höhere Näherung einer von Mott-Smith angegebenen Behandlung betrachtet 
werden. 


A. Einleitung 


Bei der Behandlung gasdynamischer Probleme mit Hilfe der Boltzman we 
schen Fundamentalgleichung?-®) ergeben sich Schwierigkeiten, wenn starke 
Abweichungen der Geschwindigkeitsverteilung der Gase von der Maxwell- 
Boltzmannschen Gleichgewichtsverteilung zu erwarten sind. In diesen 


') Vorgetragen auf der Tagung des Verbandes Deutscher Physikalischer Gesell- 
schaften 1959 in Berlin. 

2) J.H. Jeans, Dynamische Theorie der Gase. F. Vieweg Sohn, Braunschweig 
1926. 

3) A.Eucken u. K.L. Wolfs Hand- und Jahrbuch der Chem. Physik, Bd. III, 
Abschn. IV, Freie Weglänge und Transporterscheinungen in Gasen von K. Herzfeld, 
Leipzig 1939. 

4) K.Oswatitsch, Gasdynamik. Springer, Berlin 1952. 

5) A.Sommerfeld, Vorlesungen über Theoretische Physik, Bd. V, Thermodynamik 
und Statistik. Dieterichsche Verl.-Buchhandlg., Wiesbaden 1952. 

6) S.Chapman u. T.E.Cowling, The Mathematical Theory of Non-uniform 
Gases. At The University Press, Cambridge, Repr. 1953. 

7) H. G. Reik, Reibungstensor, Diffusions- und Wärmestrom in stark inhomogenen 
Gasen mit Maxwellscher Wechselwirkung. Habilitationsschrift, Aachen 1957. 

8) Handbuch der Physik. Bd. XII, Thermodynamik der Gase. Springer, Berlin-Göt- 
1958. 
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Fällen wurde bisher versucht, auf der Basis formaler Reihenansätze für die 
Geschwindigkeitsverteilung Näherungslösungen der Boltzmann-Gleichung 
zu gewinnen. Diese Ansätze sind Entwicklungen einer Maxwell-Verteilung 
nach den Komponenten der Teilchengeschwindigkeit, die der Form nad 
Hermiteschen Polynomen in drei Veränderlichen ähnlich sind, z. B, 
die Entwicklung von Burnett?)!%). Die Handhabung des Burnettschen 
Ansatzes in Verbindung mit der Boltzmann-Gleichung und im Zusammen. 
hang mit praktischen Näherungsverfahren für die Lösung der Diff.-Gleichung 
ist sehr aufwendig. Man ist dann gezwungen, sich auf wenige Naherungs. 
schritte zu beschränken und man kann nur Gasbewegungen mit relativ kleinen 
Abweichungen vom thermodynamischen Gleichgewicht behandeln. 
ersten Entwicklungsgliedern ist der Burnettsche Ansatz zu wenig flexibel, 


vorzeitig abgebrochen wird. 


Der stationäre Verdichtungsstoß a 


Die Problematik des Burnettschen Ansatzes bei stärkeren Abweichungen 
vom thermodynamischen Gleichgewicht wurde deutlich am Schulbeispiel des 
stationären Verdichtungsstoßes!!) 12), bei der Berechnung der Stoßfrontstruktur 


und Taylor (1910) zurück. Sehr viel weiter gehende Untersuchungen sind 


In den fr 


wenn stärkere Abweichungen vorliegen. Es können im übrigen auch negative, 
unphysikalische Anteileder Verteilung vorkommen, wenn die Reihenentwicklung 


und bei der Ermittlung der sich in dieser Front einstellenden Geschwindig. | 
keitsverteilung. Die ersten Stoßfrontbehandlungen gehen auf Lord Rayleigh 


F. Wei 
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von Becker!) angestellt worden. Die aus den Navier-Stokes-Gleichungen 
bei Verwendung konstanter Reibungs- und Wärmeleitungskoeffizienten ge- 
wonnenen Ergebnisse wurden später von Thomas!) verbessert. Diese Er- 
gebnisse — unter Zugrundelegung des Burnettschen Ansatzes — bezeichnet 
man mit Lösungen zweiter Näherung. Lösungen dritter Näherung sind 1951 
von Zoller?) untersucht worden. Fehlerabschätzungen zeigen dabei, daß die 
Näherungen schon bei Stoßstärken von etwa M = 2,3 nicht mehr ausreichend 
sind. Es wird auch erkennbar, daß ein nächster Näherungsschritt zu einem 
kaum noch tragbaren Rechenaufwand führt. 

Ebenfalls im Jahre 1951 wurde von Mott-Smith**) ein neuer Weg ein- 
geschlagen dadurch, daß schon bei der Wahl des Ansatzes für die Geschwindig- 
keitsverteilung apriori charakteristische Züge des vorliegenden gasdynamischen 


Math. Soc. London 39, 385—430 (1935). 

10) D. Burnett, The Distribution of Molecular Velocities and the Mean Motion in a 
non-uniform Gas. Proc. Math. Soc. London (Sec. Series) 40, 382—435 (1936). 

11) G. G. Stokes, On a Difficulty in the Theory of Sound. Philis. Mag. 33, 349—356 
(1848); Collected Papers, Band 1. 

12) E.F. Greene u. J.P. Toennies, Chemische Reaktionen in Stoßwellen. Dr. 
Dietrich Steinkopf, Darmstadt 1959 (Reihe: Fortschritte der physikalischen Chemie, 
Bd. 3). 

13) R.Becker, Stoßwelle und Detonation. Z. Physik 8, 321—362 (1923). 


Chem. Phys. 12, 449-453 (1944). 


Physic. Rev. 82, 885—892 (1951). 


®) D. Burnett, The Distribution of Velocities in a slightly non-uniform Gas. Proce. |: 


14) L.H. Thomas, Note on Becker’s Theory of the Shock Front. The Journ. of 


15) K. Zoller, Zur Struktur des Verdichtungsstoßes. Z. Physik 130, 1—38 (1951). 
16) H.M. Mott-Smith, The Solution of the Boltzmann Equation for a Shock Wave. 
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‚IWeglänge von beiden Seiten aus der Fall ist. Als Ansatz für die Geschwindig- 
"[Temperatur, d.h. mit den Werten, die das Gas vor und hinter der StoB- 


‘liberdies auf das spezielle Problem des Verdichtungsstoßes zugeschnitten und 


R 


en 


F. Weitzsch: bei vom Gleichgewicht ‚405 


Problems berücksichtigt werden. Mott-Smith geht von der Vorstellung aus, 
daß die Geschwindigkeitsverteilungen der stationären Punkte auf der Vorder- 
seite und auf der Rückseite der Stoßfront in die Stoßfront hineinreichen und 
im Stoßfrontbereich mit verschieden großer Teilchenzahl nebeneinander be- 
stehen. Man kann vermuten, daß das erstere wenigstens etwa über eine freie 


keitsverteilung wird daher die Summe zweier Maxwell-Verteilungen gewählt 
- mit jeweils konstanten Werten :der Strömungsgeschwindigkeit und der 


front besitzt. Variabel sind nur die beiden Teilchendichten. Die Lösungen 
lassen jedoch keine befriedigenden Fehlerabschätzungen zu. Der Ansatz ist 


ist daher nicht an die Stelle allgemeiner Entwicklungen (z. B. Burnettscher 
Ansatz) zu setzen. 


ßB. Formulierung eines neuen Ansatzes für die Gewinnung von Näherungs- 
lösungen der Boltzmann-Gleichung 


Ein Ansatz, der gegenüber den Reihenansätzen wesentliche Vorteile bietet, 
ist wie folgt aufgebaut. Wir denken uns zunächst eine beliebig kleine Gruppe 
von Teilchen des Gases, die eine Maxwell-Boltzmannsche En richts- 


verteilung 
m 3/2 m 


ten ge- 


ese Er- find eine mittlere Temperatur T. Wir summieren weiter über unendlich viele 
eichnet folcher Gruppen mit verschiedenen Werten » und 7 mit Hilfe einer örtlich 
id 1951 find zeitlich variablen Gewichtsfunktion g(v, T,r,t). Dabei ist zu beachten, 
daß die Haß » und 7 für die Gesamtheit ihre Bedeutung als physikalische Strömungs- 


eichend 


1 einem frariable darstellen. Damit entsteht der neue Ansatz für die Verteilungsfunk- 


Jeg ein- 
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'k Wave. 


iat. Diese Gruppe hat also eine lokale mittlere Strömungsgeschwindigkeit v 


schwindigkeit und Temperatur verlieren und daher lediglich Integrations- 


dus Normierungsgründen verlangen wir 


T>0; T+T. 
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Schließlich ist 


k 


2 
p? + 3 = 
Mit Verwendung von Gl. (3) erhält man 


=f eFde. 


T = f [7 + (o — g dv aT. 


Einige Eigenschaften des Ansatzes sind die folgenden: 
a) Es kann jede beliebige Verteilungsfunktion durch ihn dargestellt 
werden. Die Gewichtsfunktion gestattet in einer Vielzahl von Möglichkeiten 
eine rasche Anpassung der Verteilungsfunktion an eine gegebene Aufgabe. 
b) Beim Übergang von g zu einer einzelnen Diracschen 6-Funktion an den 
Stellen 


T>-T 


mit dem Inhalt 1 erhält man die Maxwell-Boltzmannsche Gleichgewichts 
verteilung. 
c) Falls g im Raum 
—co< p< +00; 0<T<co 


überall positiv ist, ist auch F überall positiv. Dies ist eine hinreichende, aber 
ae nicht notwendige Bedingung. F kann auch dann noch überall positiv sein, 

wenn es Teilgebiete von v und 7 gibt, in denen g < 0 ist. Im letzteren Falk 
muß die Prüfung, ob F überall positiv ist, nach Lösung der Bewegungsaufgabe 
erfolgen. 

d) Durch eine Summe aus n 6-Funktionen für g erhält man eine Verteilung 
aus n verschieden stark vertretenen gewöhnlichen Gleichgewichtsvertelungel 
mit voneinander verschiedenen Parametern v,, 7,. Im Fall n = 2 und wenn 
die v;, 7, konstante Größen sind, liegt der Ansatz von Mott-Smith vor 


€. Herleitung der gaskinetischen Differentialgleichungen 


Der vorliegende Ansatz ermöglicht keine direkte Integration der Boltz- 
mann-Gleichung. Da es jedoch bei den meisten physikalischen Problemen 
nur auf die Be ‚rechnung der makrophysikalischen Größen ankommt, genügt 
die sukzessive Berechnung der Momente der Verteilungsfunktion 


Q=/QFdc 


worin Q ein Polynom in den Geschwindigkeitskomponenten £, 7, € vom Grade 
n ist. Diese Momente genügen den aus der Maxwellschen Transportglei 
chung gewonnenen Differentialgleichungen, die im Falle eindimensionaler 
ebener Strömungsvorgänge ohne äußeres Kraftfeld, auf die wir uns hier 


beschränken wollen, lauten 
et (v Q) + +2 (v E Q) = = Ist 22 (6) 


mit den Stoßintegralen 
Is =» [ Q 91 Go (101 fo2 — for fos) €, p dp dy de, de, dv, dT, 

(bezügl. p, y,c, vgl. Abb. 6 im Anhang). 
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Bei einer Beschränkung auf endlich viele Momente ist das Verfahren 
zur Berechnung der Verteilungsfunktion F, bzw. der Gewichtsfunktion g, Br; 
nicht eindeutig. Es gibt eine Mannigfaltigkeit von Gewichtsfunktionen, die 
in begrenzter Näherung zu gleichen 
Lösungen führen. Ist die Zahl der noch | 
(4) mit tragbarem Aufwand zu berechnen- g 
den Momente vorgegeben, z. B. die 6 
nicht verschwindenden Momente, die 
argestelif sich für eindimensionale Gasströmungen 
ichkeitenf| bei der Verwendung von (5) bis zum 
ifgabe, | dritten Grade ergeben, dann können 
mn andenf ebensoviele Parameter der Gewichts- 

| funktion eingeführt werden. Dabei ist v 
wegen der Vieldeutigkeit der Integral- Abb. 1. Beispiel für eine 6-parametrige 
zewichts-| funktion nicht sehr kritisch. Abb. 1 


i zeigt ein Beispiel. Die Gewichtfunktion Funktionen des Ortes und der Zeit. 


besteht aus zwei Diracschen 6-Funk- 
tionen, deren Inhalte und Koordinaten in der v—7-Ebene die 6 Parameter 
darstellen und Funktionen des Ortes und der Zeit sind. 

Ein mit Verwendung der Verteilung (1) angelegtes Verfahren bietet einige 
wesentliche Vorteile. Die Stoßintegrale können allgemein berechnet werden, 
ohne daß die Gewichtsfunktion schon bestimmt werden müßte. Als Beispiel 
geben wir die zu einer 6-parametrigen Gewichtsfunktion gehörigen Gleichungen 
an. Das Gleichungssystem sollte Lösungen ergeben, die üblicherweise mit 


nde, aberf 
itiv sein, 
ren Falle 
rsaufgabe 


zZ „Lösungen dritter Näherung‘‘ bezeichnet werden. Die angegebenen Stob- 
a a integrale gelten für das Stoßmodell mit „Maxwell-Molekülen“. Der Rech- 
a nungsgang ist im Anhang 1. angedeutet. 
r Boltz- PS vg do) 
roblemaf Q = &: a8 
nil ar. kT 
m Grade al” —9.de) + | vgdo +3» | dg) =H, (10) 
sportglei- 


uns hier de) + de) = (11) 
sowie fiir Q Ei ck mit 1 +) + k 3 
(» | vg do + 3y | do) 


‘ 
é wre /k T\2 
4 Oy 3 = a 
ft = g do + g do H, 
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+22): De 
den ausgewerteten Stoßintegralen 
H, = —H, =— v8 
7 BE, ist gesetzt dvdT = do und 1,37//2m@ = J,Y2m@ = s [vgl. Gl. (37) 
ne im Anhang]. @ ist die Kraftkonstante der Maxwell-Moleküle. 


Pie a" q In den Gln. (8) bis (13) kénnen die ersten drei der sechs Parameter der 
en _ Gewichtsfunktion mit Hilfe der Integrale (2), (3) und (4) durch die physika- 
Tisch n Größen », v und 7 ersetzt werden 


fvgdo 


ze" man für die übrigen in den Gln. (8) bis (13) vorkommenden Integrale 
noch einige Ersetzungen ein 


vf (v2 — #2) gdo =—2F, (17) 
af = —40¥, —5¥, + 5%, 
(v—v)gdg = — B, 
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|| 
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dann lassen sich die Gleichungen (8) bis (13) i in eine bekannte Gestalt bringen. 
Die Y, verschwinden alle bei VorhandenSein einer Maxwell-Boltzmann- 

schen Gleichgewichtsverteilung und müssen in den Bewegungsgleichungen auch 
als Ausdrücke für die Beschreibung der Transportphänomene zu deuten sein. 
Aus den Gln. (8) bis (1 u nn man exakt 


af kP\ a 
+2 (> 2 (4m) (19) 
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Y, und Y, sind proportional den Burnett-Termen B,, und By, (vgl. °), 19) 
und 15). Verfährt man mit den übrigen drei Diff.-Gln. so, wie bei dem itera- 
tiven Lösungsverfahren mit dem Burnettschen Ansatz, in dem man an- 
nimmt, daß Y, und Y, sich wie die Burnett-Koeffizienten B,, und By in 
Reihen von wachsender Ordnung in den Ableitungen nach v und T schreiben 
lassen, dann erhält man 


- | ö | 


Man erhält mit den angeschriebenen ersten Gliedern von Y, und YW, die 
Navier-Stokes-Gleichungen. Die Koeffizienten «u und A sind die gleichen, 
wie man sie bei Verwendung des Burnettschen Ansatzes (für Maxwell- 
Moleküle) gewinnt. Dies bedeutet hier, daß sich bei dem Grade der vorliegen- 
den Näherung die Navier-Stokes-Gleichungen herleiten lassen, ohne die 
Gewichtsfunktion g bestimmen zu müssen. Daraus folgt ferner, daß beliebige 
sechs-parametrige Verteilungen zu gleichen makrophysikalischen Lösungen 
führen, sofern — vom Grade der Näherung festgelegt — die Abweichungen ag 
der Verteilung vom Gleichgewichtstyp hinreichend klein sind. tam 
Andererseits erkennt man, daß für stärkere Abweichungen die Glieder j 
höherer Ordnung von Y, und Y, sicher nicht unberücksichtigt bleiben dürfen, 
d. h., die Integrale höheren Grades in (17) müssen in der Rechnung mitge- 
führt werden. Für das praktische Rechnen ist es bequemer, nicht von den 
Diff.-Gln. (18) bis (21) auszugehen, sondern unmittelbar von dem Gleichungs- 
system (8) bis (13). Während es beim Burnettschen Ansatz naheliegt, ein a 


durch die Reihenentwicklung bedingtes iteratives Lösungsverfahren zu ver- | 
wenden, liegt hier kein Anlaß vor, die Differentialgleichungen in ein iteratives Bu: 
System näherungsweise zu verwandeln. 
D. Anwendung des neuen Ansatzes am Beispiel des stationären Verdichtungs- ia 
stoBes 
Normierungen 


Fir die Kennzeichnung der Starke eines VerdichtungsstoBes werden im 
allgemeinen die Machzahl, das Druckverhältnis oder das Dichteverhältnis 
verwendet. Besonders angenehm erweist sich im Zusammenhang mit zweck- 
mäßigen Normierungen ein Faktor 


3 


das Dichteverhältnis bedeutet. 
Der Wertebereich von rn ist 0 <¢ < 3/5; fiir 1 En <0. 
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5 
(%, 7, = Größen vor der Stoßfront) 
An den stationären Punkten vor und hinter der Stoßfront gilt dann 
Für die Ortskoordinate x führen wir ein oa 29 » 


- 


Ay = freie Weglänge vor der Stoßfront. (24) 
0 


Lésung von Mott-Smith 
Wir gehen zunächst von dem Modell von Mott-Smith!®) aus, dessen Ansatz 
im Bilde der Gewichtsfunktion die Gestalt hat 
g(u, = d(u + 1) + 1) + a,(Z) — 1) 6(8— 1). (26) 

(Mit dem Index 0 ist der Zustand vor der Stoßfront, mit dem Index 1 der 
Zustand hinter der Stoßfront bezeichnet.) Es gibt zwei freie Parameter ay und 
a,. Die Randbedingungen sind 

a)(— co) = 1 a,(— co) = 0 

(+00) = 
Mit Hilfe der Momente bis zum 2. Grade erhält man die Lösung von Mott- 
Smith 


(26) 


a= u2); u = tanh (x 2). (27) 

Die aus der Wendetangente folgende StoBfrontbreite / ist int ahs 
l= 2A/x (28) 


hängt nur von der StoBstarke ab. Es gilt fiir Ma x well-Molekiile als Funktion 
q 


‘16 107 ~ 
Fiir das Modell starr-elastischer Kugeln folgt (in einer Reihe geschrieben) 


n!(2n— 


(30) 
h of! 7 126 1386 


j nr Für die Strömungsgeschwindigkeit v und für die Temperatur 7 führen wir I 
> . . . 
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Via 5(v — 1)? 


0<R < 3/2. 


Die reziproken Stoßfrontbreiten (auf die freie Weglänge A) bezogen) sind in 
Abb. 2 als Funktion des Dichteverhältnisses 5 dargestellt. Weiterhin sind die — 
mit Hilfe des Burnettschen 


Ansatzes für Maxwell-Mole- | 
kiile sich ergebenden Lösungen Becker - Thomas | | 
zweiter Näherung (Becker- 45 
Thomas) und die drei von ) 
Zoller berechneten Werte 04 
Der letzte Wert von Zoller 03 Hs a 
bei einem Druckverhältnis Mott- mh | 
p= 6,5 bzw. 9 = 2,571 ist, 02 
wie wir bereits eingangs be- 
merkten, nicht mehr zuver- 01 
lässig. Zoller 
0 
Lösung dritter Näherung 15 20 25 30 35 9 40 


£ r ke Abb. 2. Reziproke, auf die freie Weglänge des un- 
Wir nehmen nun im näch- gestörten Gases bezogene Stoßfrontbreite als Funk- 
sten Näherungsschritt eine tion des Dichteverhältnisses (für c,/c, = 5/3). Es 
Erweiterung des Ansatzes von die 
F . . omas und dritter Näherung von Zoller (Max- 
Mott-Smith Var; in dem wit  well-Molekiile) sowie die Lösungen von Mott- 
zu einer 6-parametrigen Ge- Smith für Maxwell-Moleküle und für starrela- 
wichtsfunktion übergehen. Es stische Kugeln eingetragen. 


bietet sich an, in der Gewichts- 
funktion. (25) außer a,(2) und a,(z) auch noch die Lagekoordinaten der ö- 
Funktionen in der v—7'-Ebene bzw. u — 9-Ebene als Parameter einzuführen 


g(u, &) = a,(&) ölu — — (z)) 


(31) 
+ a, (%) ö(u — (x)) — 9, (z)) 
Die Randbedingungen sind 
Ay(— co) = 1; (— co) = U = —1; =) = —1 


a,(— co) = 0 a 
ay (+ 00) = 0 (32) 
a,(+ co) = 1; u,(+ 00) = = +1; 00) =%),= +1 


(vgl. Abb. 1). Mit Hilfe der 6 freien Parameter lassen sich die Momente bis 
zum dritten Grade, d.h. die Gln. (8) bis (13) erfüllen. 

Wir verzichten auf die Wiedergabe der Rechnung, in der keine prinzipiellen 
auftreten. Mit Hilfe der unmittelbar zu Glei- 


(23) 
ınn | 
(24) 
Ansatz 
(25) 
c 1 der 
und 
(26) 
Mott- 
4 
| 
| 
1 
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chungen wird man auf ein System von drei gewöhnlichen, nichtlinearen Diffe- 
rentialgleichungen geführt. 


Da nicht nur die makrophysikalischen Größen », ö und 7 berechnet werden 
können, sondern auch die Geschwindigkeitsverteilung selbst (im Rahmen der 
zugrunde liegenden Näherung), kann man durch Vergleich der Verteilungen 
verschiedener Lösungen die Güte der Lösungen beurteilen. 


Modell mit 
bew.- Fkt. 


2,/0, 
Abb. 3. £&-Komponente der Verteilungsfunktion verschiedener Lösungen für die Stoß- 


frontmitte bei einem starken Stoß. Die mit dem neuen Modell gewonnene Lösung (Modell 
mit Gew.-Fkt.) ergibt die glatteste Verteilung. 


Bei dem Ansatz von Mott-Smith sind die Orte der beiden 6-Funktionen 
in der v—T-Ebene bzw. u—d-Ebene per definitionem an den stationären 
Punkten unveränderlich. In unserer dritten Näherung mit variablen Koor- 

dinaten ergibt sich, daß sich 
1 die 6-Funktionen bemerkens- 
wert wenig von den statio- 


04 Mott- Smith La nären Punkten entfernen, so 
daß der bimodale Charakter 


| 


0,3 = der Verteilung erhalten bleibt. 

Aber auch schon bei Verwen- 

02 Modell dung des Burnettschen An- 
| 


satzes ist diese bimodale 
Eigenschaft sehr deutlich er- 
~ . 
2oller | kennbar. In Abb. 3 haben wir 
0 | N die &-Komponente der Ver- 
7 15 20 23 3.0 35 9 40  teilungsfunktion (F') verschie- 
Abb. 4. Reziproke, auf die freie Weglänge des un- dener Lösungen für die Mitte 
gestörten Gases bezogene StoBfrontbreite wie Abb.2. der Stoßfront konstruiert. Es 
Die für drei starke Stöße mit dem neuen Modell . t k 5 laß di it 
berechneten Werte liegen auf einem als Fortsetzung 
der Kurve durch die Zollerschen Werte dritter unserem Ansatz erhaltene Ver- 
Näherung sich zu denkenden Kurvenstückes. teilungsfunktion (Modell mit 
Gew.-Fkt.) die glatteste Ver- 
teilung liefert, obwohl sie aus zwei Verteilungsfunktionen zusammengesetzt 
ist, deren Maxima nur wenig von den stationären Punkten entfernt sind. Die 
eingangs erörterte Vorstellung von Mott-Smith erhält damit eine offensicht- 
liche Rechtfertigung. 


mit bew.- Funktion 


4 
| 
er 10 N N 


Der Verlauf der mittleren Strömungsgeschwindigkeit als Funktion des Ortes | 
ist bei den einzelnen Lösungen nur wenig voneinander verschieden. Die rezi- 
proken Stoßfrontbreiten sind in Abb. 4 aufgetragen (vgl. auch Abb. 2). Wir 
haben — da sich verhältnismäßig leicht zeigen läßt, daß bei schwachen Stößen =| 
nur kleine Abweichungen von der Lösung (27) vorliegen — drei Fälle stärkerer 


Stöße berechnet, und zwar für die Stoßstärken. oie 

(Der letzte Wert ist nur theoretisch von Bedeutung.) Die drei von uns berech- . 


neten Werte liegen in der Nähe eines als Fortsetzung der Kurve durch die 
Zollerschen Werte sich zu denkenden Kurvenstückes. ” 


10 
Abb. 5. &-Komponente der mit Hilfe des neuen Modells berechneten Verteilungsfunktion 
an verschiedenen Stellen der Stoßfront. Die Verteilungsfunktion vor der Stoßfront 
(@ = — 1) geht beim Durchtritt des Gases durch die Stoßfront (@ = Oentspricht der Mitte __ 
der Stoßfront) über in die Verteilungsfunktion auf der Rückseite der StoBfront (@ = + 1). | 


Interessant ist noch das Verhalten der Verteilungsfunktion an verschiedenen 
Stellen der Stoßfront, d. h., die Art und Weise, in der das Gas beim Durchtritt 
des Gases durch die Stoßfront seine Geschwindigkeitsverteilung ändert (Abb. 5, — 
vgl. Abb. 3). Vor der Stoßfront gilt die Kurve für @ = — 1, auf der Rückseite _ 
der StoBfront die Kurve für wu = + 1. . 


Schlußbemerkungen 

Für die Lösung gasdynamischer Probleme bei starken Abweichungen vom 
thermodynamischen Gleichgewicht haben wir einen Ansatz für die Verteilungs- 
funktion vorgeschlagen, mit Hilfe dessen man schon in den ersten Näherungs- 
schritten eine bessere Anpassung an das jeweils vorliegende physikalische _ 
Problem erreichen kann als bei der Verwendung formaler Reihenansätze. 

Am Beispiel des stationären Verdichtungsstoßes haben wir den neuen Ansatz 
erprobt. Die von uns ermittelte Stoßfrontbreite eines theoretisch unendlich Ü ER 
starken Stoßes beträgt drei freie Weglängen des ungestörten Gases vor der Er 
Stoßfront oder (mit einem Dichteverhältnis » = 4,0) zwölf freie Weglängen des _ 
Gases auf der Rückseite der Stoßfront. 


q F. Weitzsch: Behandlung gasdynamiseher Probleme bei Abweichungen vom Gleichgewicht 413 _ ai 
ffe- 
den 
der 
3 
IT - ZN 
ch 
1S- 
i0- 
so 
er 
rt. 
n- 
n- 
‚le 
T- 
ir 
T- 
e- 
be 
it 
T- 
it 
T- 
zt 
ie 
t- 
| 


414 Annalen der Physik. 7. Folge. Band 7. said etic eat 


Wir denken, daß der sehr allgemeine Ansatz fiir die Verteilung unter Ver- 


wendung einer Gewichtsfunktion sich auch auf andere Aufgaben der klassischen 
kinetischen Gastheorie anwenden läßt. 


Herrn Prof. Dr. G. Burkhardt, in dessen Institut der vorliegende Bei- 
trag entstand, sei an dieser Stelle für die vielen Anregungen und für das freund- 
liche Lehrer-Schüler-Verhältnis gedankt. Auch Herrn Dr. R. W. Larenz 
gebührt ein besonderer Dank für zahlreiche anregende Diskussionen. 


. Auswertung der Stoßintegrale 


Unter Maxwell-Molekülen versteht man Teilchen, die sich mit einer der 
fünften negativen Potenz der Entfernung proportionalen Kraft abstoßen 


K = G m? 1-5, (33) 
Einführung einer Größe (vgl. ?)) 
Cc, 
(c, = Relativgeschwindigkeit) 


ergibt sich für den Stoßwinkel + 


¢ 

a a nach Abb. 6 ein elliptisches Integral 

p 

al2 

Fr dd 
y Vı — sin? (3) 
0 

Abb. 6. Hilfsskizze zur Auswertung der mit ie 
Stoßintegrale. 1 1 re 
t= m) 


Bei der Auswertung der StoBintegrale kommt bei den Momenten bis zum 
dritten Grade der Winkel y nur in einem einzigen Zn Integral vor") 


Das Stoßintegral (7) lautet nach Einführung von « und mit einigen bekannten 


Umformungen (vgl. 5), S. 286) 
v2 wad 4 
54 = 3 v2 mG Sa 92 for fo2(@i + — Qı — (36) 
ER . de, de, dy x dx dv, dT, dv, dT». 


Führt man eine Größe 


s=J,/2mG =1,37/2mG 


(37) 


17) Das Integral wurde schon von Maxwell, Aichiund Tanakadate, sowie von 
Chapmann berechnet und von uns, um Irrtiimer durch Ersetzungen zu vermeiden, 
verifiziert (vgl. auch 2), S. 305 und "), S. 9). 


F.W 


sowi 


ein, 


ode 


win 


mit 


ein 


= 
Da 
: Ma: 
Z. B 
| Cp 
\ / 
\ 
| 
| 
Mit 
ver 
erh: 
{ sin®y (35) 
= 
= 
= 


ein, dann lautet das Stoßintegral (7) 


Da die Gewichtsfunktion g nicht von c abhängt, sind nur Integrationen über 
Maxwell-Verteilungen durchzuführen, die keine Schwierigkeiten bereiten. 
Z. B. erhält man für Q = & bei eindimensionaler Gasströmung 


—(fvgdvat)} 
oder mit Verwendung von Gl. (3) 


I«(Q= &) = — (v— dvdT. 


Starrelastische Kugeln 


Im Modell starrelastischer Kugeln kann man mit Verwendung des Stoß- 
winkels x schreiben 


c, pdpdy = ¢, 7 sin x dy dy. (41) 
Mit dem Wirkungsquerschnitt o = 42 r} (r, = Molekiilradius). Die Rechnung 


verläuft in gleicher Weise wie beim Modell mit Maxwell-Molekülen. Man 
erhält z. B. für Q = & 


32 7/k/T,+T 
Is, = — ros f 9192 )) 
5 3 = 
-dv, dT, dv, dT, 


(v, — v2)? 


k 
2) 
2. (7, + Ty) 


Der Ausdruck in der geschweiften Klammer läßt sich übersichtlicher in 
einer Reihe schreiben 

2 U" 15 

n=0 Rn! (2n —1)(2n + 1) (2n+3) (2n+5) 


4 3 
2 
Cia : F. Weitzsch: Behandlung gasdynamischer Probleme bei Abweichungen vom Gleichgewicht 415 5 val 
Ver- | sowie ein weiteres abspaltbares Integral = 
Cg) (Q Q, — Qe) ¢, p dp : 
Bei- 
und- 
St y“ dC, aly D, 1 ( Do 9° 
| 4 
| | 
der om 
n 
(33) 
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oral 
(34) 
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r!?) 
(35) 
iten = 
; P(x) 
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oder ausgeschrieben 
U + U F 


er Im Falle einer aus zwei Diracschen 6-Funktionen nach Gl. (31) bzw. Abb. | 
bestehenden Verteilungsfunktion ist die Integration der Integrale (40) und 
(42) besonders einfach. Es folgt für Maxwell-Moleküle 


Ts¢(Q = &) = — 8 ay a, (v9 — (43) 
worin jetzt die Indizes 0 und 1 die beiden ö-Funktionen bezeichnen. 
Für starrelastische Kugeln gilt : 4 
32 (To + Ts) 
(Q = #) Ay a, (Vy — 
(44) 
(Up — 
(0< R < 3/2). 
2— (To + Ti) 


In gleicher Weise lassen sich die übrigen StoBintegrale auswerten. er z 
‘ 


2. Modellvergleich weer 


Beim Modell starrelastischer Kugeln enthalten die StoBintegrale den ver- 
haltnismaBig leicht bestimmbaren Wirkungsquerschnitt o, der im Falle des 
StoBfrontproblems bei Normierung der Ortskoordinate auf die freie Weg- 
länge überdies herausfällt. Bei Maxwell-Molekülen ist die Definition eines 
Wirkungsquerschnittes bei der Ablenkung der Teilchen im Kraftfeld’ mit der 
Konstanten @ nach Gl. (33) bzw. mit dem Faktor s nach Gl. (37) ohne Willkür 
schwierig. Man ist auf abschätzende Vergleiche angewiesen. 

Ein offenbar brauchbares Verfahren besteht darin, daß die StoBintegrale 
miteinander verglichen werden, insbesondere in der Nähe der Gleichgewichts- 
zustände. Dies ist in den Gln. (43) und (44) für R<1 der Fall, so daß man un- 
mittelbar als Vergleichsformel erhält 


|Vergl. 15/2 m 2 


(45) 
T, und T, sind beim Ansatz von Mott-Smith die Temperaturen in den 
stationären Punkten, bei anderen Modellen wird man (7, + 7,)/2 als Mittel- 
wert der Temperatur ansehen können. Von Jeans (vgl. ?), S. 358) wird auf 
einem anderen — etwas schwierig zu übersehendem — Wege ein Faktor ge- 
funden ar 
8 | _ 2/2 - 1,37 ; (46) 
0 |Jeans m 
Diesen Wert hat auch Zoller von Jeans übernommen (vgl. ®), S. 31). 
Für den Fall 7 = (7, + T,)/2 erhält man mit Gl. (45) und (46) 
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und es ist zu sehen, daß beide Ergebnisse für diesen Fall als praktisch gleich 
angesehen werden können. 

Als Probe haben wir auch noch das Stoßintegral für @ = & untersucht, 
wobei sich ebenfalls die Vergleichsformel (45) ergibt. Verwendet man die Ver- 
gleichsformel (45) in der Lösung von Mott-Smith (27), dann ergibt sich für 
die Stoßfrontbreiten 

ho Ag | 1 

[vgl. die geschweifte Klammer in Gl. ia. 
worin im Extremfall des unendlich starken Stoßes — dort ist R = 3/2 — der 
Wert in der geschweiften Klammer 1,20, die maximale Abweichung der rezi- 
proken Stoßfrontbreiten also 20% beträgt (vgl. auch Abb. 2). Bei Verwendung 
des Momentes Q = & an Stelle von Q = & beträgt die maximale Abweichung 
nur 8%. In unserer Lösung mit variablen Koordinaten der 6-Funktionen 
sind die Abweichungen noch geringer. 


Hannover, Institut für Theoretische Physik der Technischen Hochschule. 
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Über den Beitrag der Photonen zur Elektronenauslösung 
aus der Kathode in Glimmentladungen in Hg-Dämpfen 


dm: 


Von E. Badareu und F. Wachter 4 


Inhaltsübersicht 


Es wird eine Versuchsanordnung angegeben, die es erlaubt, den Beitrag des 
P’hotoeffektes einer anomalen einfachen Glimmentladung und einer Doppel- 
kathodenentladung zur Elektronenauslösung an der Kathode zu bestimmen. 
Dieser Beitrag wird für verschiedene Kathodenfälle in Hg-Dämpfen bestimmt. 

Aus den MeBergebnissen kann gefolgert werden, daß der Photoeffekt in 

_ kleinerem Maße zur Elektronenauslösung beiträgt, kaum aber zur Erklärung 
2 der Stromamplifikation, durch die der Hohlkathodeneffekt charakterisiert 
wird, werden kann. 


Einleitung 


In einer Untersuchung des Hohlkathodeneffektes, oder auch Doppel- 
kathodeneffektes, wie derselbe bei ebenen, parallelen Kathoden noch genannt 
werden kann, kommen Little und von Engel!) zur Schlußfolgerung, daß 
der Photoeffekt, dem auch bei der anomalen Glimmentladung der Haupt- 
anteil an der Elektronenbefreiung aus der Kathode zugeschrieben wird, wesent- 
lich zu der bei diesem Effekt beobachteten Stromverstärkung beiträgt. 

In einer früheren Arbeit haben wir?) zunächst den Doppelkathodeneffekt 
in Hg-Dämpfen gemessen und festgestellt, daß er große Werte erreicht. Würde 
der Photoeffekt wesentlich hierzu beitragen, so müßte sich das auch in einer 
entsprechenden Vergrößerung desselben, gegenüber der einfachen anomalen 
Glimmentladung widerspiegeln. 

Wir konnten damals feststellen, daß es in Hg-Dämpfen, bei geeignetem 
Druck und Abstand, zu einer zusätzlichen Stromverstärkung gegenüber den 
anderen Dappetinthodens ntladungen kommt. py Intensitätskurve I bei 
V = konst. und a = konst. (V = Kathodenfall, a = Abstand zwischen den 
beiden Kathoden) verläuft nicht re eier wie es in erster Näherung bei den 
Edelgasen*), H,, O, und N, der Fall ist, sondern weist einen starken Strom- 
anstieg gegen kleine Drucke hin an. 

Ziel der Untersuchung war: Messung der Elektronenauslösung durch Ent- 
ladungsphotonen und Bestimmung des Beitrages des Photoeffektes zur Elek- 


1) P.F. Little and A. von Engel, Proc. Roy. Soc. (London) A, 224, 209 (1954). 


2) EK. Badareuu.F. Wächter, Journ. of Electr. and Control, vol.4, No. 6, 539 (1958). 
3) A. Güntherschulze, Z. techn. Physik 11, 49 (1930). 
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tronenauslösung aus der Kathode an einer stationären anomalen gewöhnlichen 
Entladung und an einer Doppelkathodenentladung, sowie Feststellung eines 
eventuellen Beitrages des Effektes zur zusätzlichen Stromverstärkung bei 
Hg-Dämpfen. 

Zur Versuchsanordnung 
Die Versuchsanordnung muß so beschaffen sein, daß sie die im Fallraum 
und im negativen Glimmlicht erzeugten Photonen möglichst ungehindert zur 
Meßstelle gelangen läßt, ohne daß auch Ladungsträger dahin kommen. Wir 
haben das in Anlehnung an eine Arbeit von Costa‘) mit Hilfe eines ,,Konden- 
sators‘‘ zu erreichen versucht. 
Die Versuchsanordnung zur Erzeugung der Doppelkathodenentladung 
gleicht unserer früher angewendeten?). Zusätzlich wurde noch eine Anordnung 
zur Messung der Photoströme in die Entladungsröhre eingebaut. 


Cc Go 

Abb. 1 


Die Elektrodenanordnung zeigt Abb. 1. Die beiden ebenen, kreisférmigen 
Eisenkathoden mit einem Durchmesser von 3,6 cm stehen sich in einer Ent- 
fernung von 0,65 cm gegenüber. Um eine Entladung auf der Rückseite zu 
verhindern, wurden sie mit je einem Glastrichter überdeckt. An den seitlichen 
Wänden dieses Trichters werden je zwei ringförmige Anoden A, und A, so 
befestigt, daß sie auf dieselbe Höhe wie die beiden Kathoden zu liegen kommen. 
Wir erreichen dadurch größtmöglichste Symmetrie der Entladung. Die beiden 
Elektrodenpaare konnten dieses Mal nicht bewegt werden. Gleich an die 
Anoden schloß sich der Kondensator C an, dessen Platten von 10”? cm Dicke 
und 3,2 em Länge, in 0,2 cm Abstand so aufgeladen werden konnten, daß auch 
ein Teilchen mit größtmöglicher Energie (bis 800 eV) nicht hindurchfliegen 
konnte. Jenseits des Kondensators befand sich eine Glasröhre mit einem 
Durchmesser von 1,8cm. Abgedeckt wurde diese Glasröhre mit einem Git- 
ter G,, darauf folgte in der Röhre in einem Abstand von 0,5 cm, ein fein- 
maschiges zweites Gitter aus Nickeldraht und dann eine aus dem Kathodenma- 
terial angefertigte ,,Photoelektrode“. 

Ein Teil der in der Entladung zwischen den Kathoden erzeugten und seit- 
lich ausgesandten Photonen konnte so, nachdem sie Kondensator und Gitter 
passiert hatten, auf die ,, Photoelektrode“‘ auftreffen und von da Elektronen aus- 
lösen, die dann ihrerseits gegen @, beschleunigt wurden und mit Hilfe eines 


4) H. Costa, Z. Physik 113, 531 (1939); Z. Physik 116, 508. 
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Galvanometers (6 - 10-1! A/mm.m) gemessen werden konnten. Der Kondensator 

und die Gitter hatten die Aufgabe, eine Diffusion der zwischen den Kathoden 

erzeugten Ladungsträger zur ‚„Photoelektrode‘‘ zu verhindern. 

Ausgeheizt wurde die aus Pyrexglas bestehende Entladungsröhre im Hoch- 

vakuum 48 Stunden lang bei 420°, dann wurde das Quecksilber eingeführt. 
Das elektrische Schaltschema zeigt Abb. 2. 
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Die Entladungsröhre befand sich in einem Thermostaten, der zuerstaufhöhere 
Temperatur aufgeheizt wurde. Eine Meßreihe begann stets bei den hohen 
Temperaturen und endete bei den tiefstmöglichen. Zur Druckbestimmung 
verwendeten wir die für diese Temperaturbereiche gültigen Verfahren. 

Der so bestimmte Hg-Dampfdruck wurde nachher auf 0 °C reduziert. 

Wir haben bei so großen Drucken zu messen begonnen, daß wir es bei der 
gegebenen Entfernung zwischen den beiden Kathoden anfangs praktisch mit 
zwei voneinander unabhängigen Entladungen zu tun hatten. So konnten 
anfangs nur die Photonen der einen Seite die Entladung der anderen Seite 
beeinflussen. Einige um die Röhre von außen angebrachte, geerdete Aquadak- 
ringe beseitigten eventuelle Kriechströme. 


Messungen 

Es wurde bei konstantem V, der Entladungsstrom jeder einzelnen Kathode 
und auch der beider gleichzeitig bestimmt (Doppelkathodeneffekt); dabei 
wurde auch der entsprechende Photostrom an der Photoelektrode gemessen. 
Infolge der groBen Periode des Galvanometers muBte auf dem absteigenden 
Ast der einfachen anomalen Entladung und der Doppelentladung eine all- 
mähliche Intensitätsverkleinerung des Stromes, auf dem ansteigenden Ast der 
Doppelentladung eine entsprechende Vergrößerung (infolge der Dichtever- 
kleinerung des Hg-Dampfes) in Kauf genommen werden. Entsprechend war 
auch der Abfall bzw. Anstieg des Photonenstromes. Praktisch nicht ver- 
ändert hatte sich das Verhältnis 
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Alle diese Ströme konnten gemessen und damit auch r, und r, bestimmt 
werden. 


Bei der Bestimmung der Potentialdifferenz zwischen ‚‚Photokathode‘‘ und 
dem zweiten Gitter mußte einerseits die teilweise Rückdiffusion der Elektronen 
zur Photokathode°)®)’) berücksichtigt und andererseits eine Amplifikation 
der Elektrizitätsträger durch Elektronenstöße verhindert werden. Zu diesem 
Zwecke wurde die Potentialdifferenz V, = 12 Volt gewählt. E/p ist dadurch 
groß genug, um eine Rückdiffusion klein zu halten; andererseits haben wir 
bei dem kleinen V, xd ~ 0, so daß der Photostrom dadurch weder verkleinert 
noch vergrößert werden konnte. 


Ergebnisse 
Betrachten wir eine Glimmentladung der Stromstärke 7. PERLE N 


I 
In einer Sekunde gelangen = Elektronen aus dieser Entladung zur Anode. 


(e = Ladung des Elektrons). Wenn jedes dieser Elektronen n’ (x, y, z, v) dv 
Photonen der Energie h v im Volumelement dv = dx dy dz an der Stelle x, y, z 
des Entladungsraumes erzeugt, so haben wir fiir den durch Photonen aller 
Energien aus der Kathode ausgelösten Primärelektronenstrom J[ py: 


= € f F(x, y,z) n’(x, y,2,»v) o,(v) dv (2) 
v v 


Die Summe steht wegen der verschiedenen Frequenzen der im Entladungs- 
raum entstehenden Photonen und das Integral wegen der Ortsabhängigkeit 
des Raumwinkels. In (2) ist 


F(x, y,z) = Wahrscheinlichkeit, daß ein im Punkt x, y,z entstandenes 

Photon auf die Kathode auftrifft (bei isotroper Strahlung ist F = 2 2 = 

taumwinkel, unter dem im betreffenden Punkt die Kathode gesehen wird). 


@,(v) = Wahrscheinlichkeit, daß ein auf die Kathode auftreffendes Photon 
der Energie hy ein Elektron auslöst. 


Mit dem gemittelten Raumwinkel F, definiert durch a ri 
F N nv) ob) J F(x, y,2,») n' (x, y, 2, v) w,(v) de 
i 
n(v) = | n' (x, y, z,v) dv eo 


die Anzahl aller im Entladungsraum von einem Elektron gebildeten Photonen er 1 
der Frequenz » bedeutet, läßt sich (2) schreiben: ; 


Torn = IF a n (v) w(v) = IF W.: 


5) N.E.Bradbury, Physic. Rev. 40, 980 (1932). 
6) L.A. Joung and Bradbury, Physic. Rev. 48, 34 (1933). 
7) M.Cravath, Physic. Rev. 45, 138 (1934). 


Pau) 
s 


Für die Photokathode können wir einen (2) entsprechenden Ansatz machen: 
w(vy) = Tk fw (3) 


«(v) = Wahrscheinlichkeit, daß ein auf die ,,Photokathode* auftreffendes 
Photon ein Elektron auslöst. 

k.f = Wahrscheinlichkeit, daß ein im Entladungsraum gebildetes Photon 
bei Voraussetzung isotroper Strahlung die ,, Photokathode“ erreicht. 

k = Freie Öffnung der beiden Gitter und des Kondensators. 

f = Durch 42 dividierter Raumwinkel, unter dem die Photokathode 
vom Entladungszentrum aus erscheint. 

= Absorptionsfaktor. 

l = Entfernung zwischen dem Entladungszentrum und der Photo- 


kathode. 
N Pa 
a= 065 cn 
| V=500 Voit = be 
= re a W- 
| en 
p 2 pt 
0 04 08 12 16 20 24 FR 0% 06 2 16 20 Pı 
Abb. 3 ns Abb. 4 
de 
In Abb. 3, 4, 5 und 6ist W als Funktion des Druckes dargestellt. Gemessen | als 
wurde r, (© = 1,2) bei 500, 600, 700 und 800 Volt. W, bezieht sich jedesmal | el 
auf die gemittelte einfache anomale Entladung und W, auf die dieser ent- q 
sprechenden Doppelentladung. 
; Die W, (i = 1,2) wurden mit Hilfe von (4) und den experimentellen | *" 
_MeBwerten aus (1) berechnet. 
sic 
7 wurde aus den geometrischen Verhältnissen zu + = 500 berechnet. ka 


1 wurde wie folgt g : 
j; wurde wie folgt gemessen: E 


An Stelle der Doppelkathode wurde eine ihrer Entladungsform ahnliche wi 
W 
B 
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Wir haben dann: 


3 
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Abb. 5 
mg Deutung der Ergebnisse 


Qualitativ lassen sich aus dem Verlauf der Kurven für unsere Entladungs- 
bedingungen in Hg-Dämpfen folgende Schlußfolgerungen ziehen: 


a) der Beitrag der Photonen zu der Auslösung der Primärelektronen exi- 
stiert und ist nicht vernachlässigbar; 

b) er wächst sowohl bei der einfachen anomalen, als auch bei der Doppel- 
entladung mit der Potentialdifferenz; 

c) gegen kleine Drucke hin wächst der Beitrag der Photoelektronen zum 
Primärstrom; 

d) dieser Beitrag ist bis zu einem gewissen Druck, der den Punkten P 
der Abb. 3, 4, 5 und 6 entspricht, bei der einfachen anomalen Entladung größer 
als bei der Doppelentladung, die z. B. für V = 700 Volt in diesem Punkt schon 
eine fünffache und für V = 800 Volt eine zehnfache Stromamplikation 


5 1, + I, aufweist. Daraus folgt zunächst sogleich, daß die Photoelektronen- 
auslösung bis hierher gar nicht zu dem Doppelkathodeneffekt beiträgt, wenn 
man von der Vergrößerung (praktisch: Verdoppelung) des Raumwinkels ab- 
sieht, derzufolge eine konstante, wie später gezeigt werden soll, kleine Ampli- 
kation auftritt. 

Nur jenseits dieses Punktes ist bei der Doppelkathode ein erhöhter Beitrag 
des Photoeffektes zur Primärelektronenauslösung über denjenigen, der durch 
Raumwinkelvergrößerung entsteht, zu bemerken. Auf den Punkt P bezogen, 
wächst W bei den verschiedenen Spannungen höchstens auf den doppelten 
Wert. Dieser Anstieg reicht, die Größe des Beitrages der Photonen mit in 
Betracht gezogen, keineswegs aus, um den Anstieg von q zu erklären. 

Es interessiert nun W,, die Wahrse sheinlichkeit, daB alle von einem 
Elektron Ges Entladungsraumes im Mittel ein Elektron 
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er ei aus einem der Kathode gleichen Material (durch den Index k gekennzeichnet) 
auslösen. 
Von W aus (4) können wir auf W, schließen, wenn: 
=, a) @(v) gleich ist mit @,(y), wenn mit anderen Worten die Photokathode 
7 aus dem Kathodenmaterial angefertigt ist, was bei uns der Fall ist, und wenn 
sich die Oberflächenbeschaffenheit der Kathode in bezug auf den Photo. 


b) wenn u(y) bei nicht vernachlässigbarem &(v) oder umgekehrt (vy), 
bei nicht vernachlässigbarem u(v) vernachlässigt werden kann. Ersteres ist 
bei den Nichtresonanzlinien der Fall, letzteres könnte bei den Resonanzlinien 
eintreffen. 
Mit o,(v) = w(y) ist 
Ion = IF Wert», 

wobei u einen gemittelten Absorptionsfaktor bedeuten soll. 

Nehmen wir an, daß im Punkte P der Abb. 3, 4, 5und 6 up ist, und 
setzen das diesem Punkt entsprechende W, in (5) ein, so erhalten wir mit 

I = e*4(1 +6) Ip 


(/o = Primarelektronenstrom ; = Plasmawirkungsgrad.) 
_ fiir den durch Photonen ausgelösten Bruchteil x des Primärelektronen- 
4 stromes einer einfachen anormalen Entladung: 


= um _ + 6)W 
0 
Unter unseren Bedingungen haben wir F ~ 0,47. Ba 


Nicht bekannt ist 6. Setzen wir F(1 + 6) = 1 für die einfache Entladung, 
so folgt 


Es ergibt sich so mit x aus 8) x ~ 101. 
ii 


Bes 8) E. Badareu u. G. Bratescu, Bul. Soc. Rom. 45, 82 (1944). 
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Die Berechnung der Wechselwirkungsenergie 
zweier gleichartiger einfach negativ geladener Ionen 
mit Elektronenstruktur 


In dieser Arbeit leiten wir eine Formel fiir die Wechselwirkungsenergie 
zweier gleichartiger, einfach negativ geladener Ionen mit edelgasartiger 
Elektronenstruktur für das Fermi- Amaldische Modell ab. Wir benutzen 
hier die Jensensche Methode, welche Gombäs auf eine einfachere Gestalt 
gebracht hat. Die weitgehenden Vereinfachungen von Gombäs erlauben es uns, 
einen analytischen Ausdruck für die oben erwähnte Wechselwirkungsenergie 
im Fermi-Amaldischen Modell abzuleiten. 


In der ursprünglichen statistischen Theorie von Thomas und Fermi 
können, wie bekannt, einfach negative Ionen nicht stabil sein. Man kann diesen 
Mangel des Thomas-Fermischen Modells nach Fermi und Amaldi durch 
eine elementare Korrektion beheben. Bezüglich der Fermi-Amaldischen 
Korrektion kann man folgendes feststellen!): Der Verlauf der Elektronen- 
dichte wird beim neutralen Atom im Vergleich zum ursprünglichen Thomas- 
Fermischen Modell verbessert. Durch das Abbrechen der Elektronendichte 
bei einem endlichen x,-Wert wird der Ladungsverlauf in den äußeren Gebieten 
des Atoms an den wirklichen besser angepaßt; bei positiven Ionen wird der 
Ladungsverlauf durch die Kontraktion der Elektronenwolke ebenfalls ver- 
bessert. Wesentlich?) ist jedoch für das Fermi-Amaldische Modell, daß 
auch einfach negative Ionen sich als stabil erweisen. Bei den einfach negativen 
Fermi-Amaldischen Ionen läuft die Elektronendichte bis ins Unendliche 
aus. Für die Elektronendichte o(r) in Atomen und Ionen ergibt sich im Fermi- 
Amaldischen Modell folgende Formel: 


mit 2 und u*= mls (vi) 


* 
In dieser Formel für o(r) ist 2 der in Einheiten u* gemessene Abstand vom 
Kern. Zund N bedeuten die Ordnungszahl und die Elektronenzahl der Atome 


1) P. Gombas, Die statistische Theorie des Atoms und ihre Anwendungen. Wien: 
Springer-Verlag, 1949. 
2) Man vergleiche hier '). ge 
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und Ionen. a, ist der Radius der kleinsten Wasserstoffbahn in der Bohrschen 
Theorie, und f(x) ist die Thomas-Fermische Funktion, welche folgenden 
Randbedingungen genügen muß: 


fo(9) = 1, fo(%) = und ap fo(%) = (2) 


Die Wechselwirkungsenergie u zweier Edelgasatome oder Ionen mit edelgas- 
artiger Elektronenstruktur ist nach Jensen?) durch folgende Formel ge- 
geben: 


.— 


+4,. (3) 


Die in dieser Formel auftretenden Energieterme haben folgende Bedeutung: 
u, ist die elektrostatische Coulombsche Wechselwirkungsenergie der punkt- 
förmigen Ionenladungen, die Energie u,, resultiert aus dem Eindringen des 
Kerns eines Atoms oder Ions in die Elektronenwolke des anderen Atoms oder 
Ions. Die Energieterme u, und u, geben uns die Wechselwirkungsenergie 
und die kinetische Energieänderung der Elektronenwolke zweier Atome oder 
Ionen; u, bzw. u,, bedeutet die Veränderung der Austauschenergie bzw. der 
Korrelationsenergie. Das letzte Glied w,, in (3) wurde von Gombias‘4) ein- 
geführt und dient zur Angabe der Korrelationskorrektion der Coulombschen 
Wechselwirkung zwischen den Elektronen. Die hier kurz besprochenen Energie- 
terme haben folgende mathematische Darstellung: 


„e[X, yı(ö) + N 72(6) + dv + ad) (6) 


D 

U. = x J + = 03°] dv = (7) 


die Elektronenzahl dieser Atome oder Ionen. Die Ladungsdichten 0, und 0, 
beziehen sich auf Atome und Ionen. y,(r,) 7; ist die mit der Ionenladung redu- 
zierte, in der Entfernungr, vom i-ten Kern wirksame Kernladung des - ten Kerns. 
Die Konstanten x,, x, und A, sind: 
= 2,871 e ay, x, = 0,7385e2 und A, = 0,093 (10) 
Für u, + u, läßt sich, wie Gombäs gezeigt hat, folgender analytischer Nähe- 
rungsausdruck herleiten, der für 6 > 4a, Gültigkeit hat ): 


= — [No + N, Ya(6)] + (a4 + 02 (0) + Bo 0, (9), (11) 
3) H. Jensen, Z. Physik 77, 722 (1932); W. Lenz, Z. Physik 77, 722 (1932). 

4) P.Gombäs, Z. Physik 121, 523 (1943). 
5) P.Gombäs, Z. Physik 98, 378 (1935). 
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wo a,, ß} und ß, Konstanten sind und folgende Bedeutungen haben: 


i 


a, = 47e dr (12a) 
0 
und 


202 


Die Indices 1 und 2 RE das erste und das zweite Atom oder Ion. 
Falls man von den Gliedern uw, und w,, in (3) absieht, und indem wir uns wieder- 
um auf die Wechselwirkung einfach negativer und gleichartiger Ionen beschrän- 
ken, für welche Z, — N = — 1 und Z, — Z, = — 1 sowie Z, = Z, = Z und 
N, = N, = N gilt, erhalten wir ß, = af =ß. “Für die W ar rgie 
bekommen wir in diesem Falle den Ausdruck: 


u=-u,+u, +4,+4, = 2y(d)e+ (4 + o0(0). (13) 


In dieser Formel sind x und f durch (12a) und (12b) gegeben, und y(r)r 
bezeichnet die mit der Ionenladung reduzierte, in der Entfernung r vom Kern 
wirksame Kernladung des Kernes. Im Fermi-Amaldischen Modell zeigt sich, 
daB einfach geladene negative Ionen stabil werden, und zwar verhalten sich diese 
wie im nicht-korrigierten urspriinglichen Thomas-Fermischen Modell der 
neutralen Atome. Das Gesamtpotential des einfach geladenen negativen Ions im 
Fermi-Amaldischen Modell kann mit Hilfe der gut tabellierten Fermi - 
schen Funktion des freien neutralen Atoms folgendermaBen dargestellt werden : 


V(r) = = 70". (14) 


r 


Die Formel (13) läßt sich mit Hilfe von{(12a), (12b) und (14) wie folgt darstellen 


4a 


4nZe , 40027, ZN: r 


In dieser Formel ist das auftretende Integral eine reine Zahl. Die letzte 
Formel zeigt, daß die Wechselwirkungsenergie u zweier gleichartiger einfach 
negativ geladener Ionen mit edelgasartiger Elektronenstruktur im Fermi- 
Amaldise hen Modell eine Funktion von 6-und Z ist. Die gut tabellierte 
Thomas-Fermische Funktion®) des freien neutralen Atoms erlaubt es uns, die 
Wechselwirkungsenergie « mit Hilfe der Formel (15) leicht zu berechnen. 
Der Autor dieser Arbeit hat in einer früheren Arbeit gezeigt, daß folgender 
Näherungsausdruck 

fo(x) = mit a = 0,7105 b= 0,03919 (16) 


6) S. Kobayashi u. T. Taima, Memoire of Faculty of Liberal Arts und Education. 
Kagawa, Japan (1956). 
7) T. Tietz, Tl Nuovo Cimento 4, 1192 (1% 56). 
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sehr gut die exakte Thomas-Fermische Funktion /,(z) approximiert. Mit 
Hilfe von (16) bekommen wir: 

Die zwei letzten Formeln (16) und (17) erlauben es uns, die Wechselwirkungs- 
energie u nach Gl]. (15) zu berechnen. In einer kürzlich erschienenen Arbeit 
hat Gaspär mit Hilfe der Gleichungen (4) bis (8) die Wechselwirkungsenergie 
u zweier gleichartiger freier neutraler Atome mit edelgasartiger Elektronen- 
struktur berechnet. Die eben erwähnte Arbeit von Gaspär gab dem Autor 
dieser Arbeit den Ansporn für die Berechnung der Wechselwirkungsenergie u 
zweier gleichartiger einfach negativ geladener Ionen [ GI. (15)] mit derselben 
Elektronenstruktur. 


Ich habe noch die angenehme Pflicht, Herrn Prof. Dr. P. Gombäs für 
die Möglichkeit zu danken, sein Theoretisch-Physikalisches Institut in Buda- 
pest besuchen zu können. Herrn Prof. Dr. R. Gaspar danke ich hier für seine 
wertvollen Ratschläge sowie auch für die Möglichkeit, sein Institut in Debrecen 
besuchen zu können. 


8) R. Gaspar, Acta Phys. Hung. 11, 71 (1960). hl ei 4 
Gala 


Lödz (Polen), Institut der Theoretischen Physik der Universität. 


re Bei der Redaktion eingegangen am 9. Juni 1960. 
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Herausgegeben von der 
Deutschen Akademie der Naturforscher Leopoldina 
Lieferbare Titel 
Nummer 5: 
Wissenschaftliche Selbstbiographie 
Von MAX PLANCK 
3. Auflage 
34 Seiten mit einem Bildnis. 1955. 8°. DM 1,35 
Orion: Über siebzig Jahre Geschichte der Physik werden lebendig, und inmitten dieser 
Zeitspanne die Geburtsstunde der Quantentheorie. Hinter jedem der bescheidenen Sätze 
dieses in seiner kristallenen Klarheit meisterlichen Selbstzeugnisses spürt man den großen 


Forscher. Das Heft wird für alle Zeiten eines der klassischen Dokumente der Natur- 
forschung bleiben. 


Nummer 6: 


. 
Johann Friedrich Naumann 
der Altmeister der deutschen Vogelkunde 
Sein Leben und seine Werke 


Von P. THOMSEN 

Nach des Verfassers Tode bearbeitet und ergänzt 

von Prof. Dr. E. STRESEMANN, Berlin. 
212 Seiten mit 12 Abbildungen im Text. 1957. 8°. Leinen DM 13.— 

Die Vogelwelt: Für die ältere Ornithologengeneration ist der Name Naumann noch heute 

ein Begriff. Wenn auch die wenigsten mit ihm durch die Lektiire seines originalen Werkes 

vertraut sein mögen, so kennen sie doch die kurz nach der Jahrhundertwende erschienene 

Neubearbeitung, den „Neuen Naumann‘, Die Biographie ist ein würdiges Mal für einen 
großen deutschen Ornithologen, dessen Wirken weit über seine Zeit rn 


- e 
Leben und Werk eines deutschen Biologen 
Von Prof. Dr. R. WEISSENBERG, Philadelphia 
63 Seiten mit 8 Bildnissen u. 1 Abbildung im Text. 1959. DM 3,60 
Durch die Veröffentlichung der Biographie von Oscar Hertwig soll die Erinnerung an 
sein Lebensbild und sein Lebenswerk wachgerufen werden, verdankt ihm doch eine ganze 
Generation von Biologen und Ärzten, die entweder unmittelbar seine Hörer und Schüler 
waren oder aus seinen Büchern und Schriften lernten, nicht nureine Bereicherung ihres 
Wissens, sondern auch die Formung ihres biologischen Weltbildes. 


Nummer 8; 
Friedrich Traugott Kiitzing 
1807-1893 
Aufzeichnungen und Erinnerungen 


Herausgegeben 
von R. H. WALTHER MÜLLER und RUDOLPH ZAUNICK, Halle/Saale 


300 Seiten mit 2 Bildnissen, 4 Abbildungen und einer Karte. 1960. 8°. 
Leinen DM 16,60 


Die hier vorliegenden ‚.Erinnerungen und Aufzeichnungen‘ des weiland Professors am 
Nordhäuser Realgymnasium, Friedrich Traugott Kützing, bringen Material in reicher 
Fülle und fesselnder Darstellung. Diese Selbstbiographie läßt uns den Lebenslauf eines 
Menschen miterleben, der, aus der kinderreichen Familie eines Müllers in Ritteburg an 
der Unstrut stammend, sich schon als blutjunger Apothekerlehrling und -gehilfe zum 
Naturforscher entwickelte, der schon als Student der Pharmazie und Naturwissenschaften 
in Halle sich die ersten Sporen als systematischer Botaniker verdiente, und der neben 
seinem Amt als Lehrer der Chemie und Naturgeschichte als Systematiker der Diatomeen 
und Algen grundlegende Werke schuf. 


JOHANN AMBROSIUS BARTH / VERLAG / LEIPZIG 
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af Meteoric Currents 

Von Prof. Dr. C. HOFFMEISTER ae 
Sternwarte der Deutschen Akademie der Wissenschaften zu Berlin in Sonneberg (Thiir.) 
VIII, 286 Seiten mit 109 Abbildungen im Text und 1 Ausklapptafel. Mit Zu- 
sammenfassungen in englischer Sprache. 1948. gr. 8°. DM 18.—, Hin. DM 20.— 


Astronomische Nachrichten: Die Darstellung folgt im wesentlichen den ver- 
schiedenen Phasen der fortschreitenden Lösung des Problems durch Sammlung von Be- 
obachtungsmaterial und dessen kritische, durch theoretische Betrachtungen geleitete Ver- 
arbeitung. — Insgesamt bietet das Buch eine große Fülle von Originalmaterial, von kri- 
tischen Diskussionen und neuen Ergebnissen. J, Wempe 


Grundtatsachen der sphärischen Astronomie 


en Von Prof. Dr. JULIUS DICK 


Wissenschaftlicher Leiter der Sternwarte Berlin- Babelsberg 
IV, 103 Seiten mit 48 Abbildungen im Text. 1956. 8°. 
DM 5.—, Halbleinen DM 6.— 


Mit diesem Buch gibt der Verfasser eine leichtfaBliche Ersteinführung in die Grundtat- 
sachen der sphärischen Astronomie und, darin einbegriffen, in die Reduktion astrono- 
mischer Beobachtungen. Stoffauswahl wie Darstellungsform sind darauf ausgerichtet, für 
die Himmelskunde Interessierte anzusprechen, sie dureh die Schärfe der Definitionen sowie 
die zwingende Logik der Deduktionen der klaseisehemr Astronomie zu fesseln und so zu 
späterem Studium tiefergehender Werke anzuregen. 

Dem Charakter des Büchleins entsprechend sind die oft in der sphärischen Astronomie 


üblichen Formelableitungen unter Heranziehung von Differentialausdrücken und Koor- 

dinatenumformungen vermieden. Es wird nur die Kenntnis trigonometrischer Funktionen 

vorausgesetzt, und die erforderlichen Ableitungen werden allein mit den Mittein der Sphärik 


Statistik und Physik der K 


Von Dr. NIKOLAUS B. RICHTER 
Sternwarte der Deutschen Akademie der Wissenschaften zu Berlin in Sonneberg (Thür.) 
VI, 142 Seiten mit 58 Abbildungen im Text. 1954. gr. 8°. DM 13.40 


Orion: Das Buch soll dem Leser einen genauen Bericht über den Stand der heutigen 
Kometenforschung geben. Die beigefügten Tabellen und Diagramme werden jedem, der 
sich mit Kometenforschung ernstlich befassen will, weiter helfen. Der Verfasser hat außer- 
dem noch einige bisher von ihm nicht veröffentlichte Ergebnisse mit in das Buch hinein- 
gearbeitet, das durch ein sehr ausführliches Literaturverzeichnis abgeschlossen wird. 


JOHANN AMBROSIUS BARTH/VERLAG/LEIPZIG 
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